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Superconductivity in cuprate superconductors occurs upon charge-carrier doping 
Mott insulators, where a central question is what mechanism causes the loss of elec¬ 
trical resistance below the superconducting transition temperature? In this review, 
we attempt to summarize the basic idea of the kinetic-energy driven superconduct¬ 
ing mechanism in the description of superconductivity in cuprate superconductors. 

The mechanism of the kinetic-energy driven superconductivity is purely electronic 
without phonons, where the charge-carrier pairing interaction in the particle-particle 
channel arises directly from the kinetic energy by the exchange of spin excitations 
in the higher powers of the doping concentration. This kinetic-energy driven d-wave 
superconducting-state is controlled by both the superconducting gap and quasipar¬ 
ticle coherence, which leads to that the maximal superconducting transition temper¬ 
ature occurs around the optimal doping, and then decreases in both the underdoped 
and overdoped regimes. In particular, the same charge-carrier interaction mediated 
by spin excitations that induces the superconducting-state in the particle-particle 
channel also generates the normal-state pseudogap state in the particle-hole chan¬ 
nel. The normal-state pseudogap crossover temperature is much larger than the 
superconducting transition temperature in the underdoped and optimally doped 
regimes, and then monotonically decreases upon the increase of doping, eventually 
disappearing together with superconductivity at the end of the superconducting 
dome. This kinetic-energy driven superconducting mechanism also indicates that 
the strong electron correlation favors superconductivity, since the main ingredient 
is identified into a charge-carrier pairing mechanism not from the external degree of 
freedom such as the phonon but rather solely from the internal spin degree of free¬ 
dom of the electron. The typical properties of cuprate superconductors discussed 
within the framework of the kinetic-energy driven superconducting mechanism are 
also reviewed. 


1. Introduction 


After intensive investigations over more than two decades, it has become clear that 
cuprate superconductors are among the most complicated systems studied in con¬ 
densed matter physics The complications arise mainly from that the parent 

compounds of cuprate superconductors are a form of non-conductor called a Mott 
insulator with an antiferromagnetic (AF) long-range order (AFLRO) [4||^, where 
a single common feature in the layered crystal structure is the presence of one 
to several Cu 02 planes in the unit cell Bi- Inelastic neutron scattering (INS) 
experiments show that the low-energy spin excitations in these parent compounds 
are well described by an AF Heisenberg model [^11] with the magnetic exchange 
coupling constant J ~ 0.1 eV. When these Cu02 planes are doped with charge 
carriers, the AFLRO phase subsides and superconductivity emerges leaving the 
AF short-range order (AFSRO) correlation still intact [m^. Although there are 
hundreds of cuprate superconducting (SC) compounds, they all fit into a universal 
phase diagram 12 as schematically illustrated in Fig. where the physical prop¬ 
erties mainly depend on the extent of doping, and the regimes have been classified 
into the underdoped, optimally doped, and overdoped, respectively. After AFLRO 
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Figure 1: (Color) Schematic phase diagram of cuprate superconductors (AFLRO region 
AF, SC-state SC, pseudogap phase PC, normal-metal phase NM, magnitude of 
the pseudogap Apc, coupling strength A*). Tn is the Neel temperature, while 
the temperature below which superconductivity (a pseudogap) is observed is 
denoted by Tc (T*). [From Ref. 
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is destroyed rapidly by doping, there are three apparent regions of the phase di¬ 
agram: (a) a d-wave SC phase, where the maximal SC transition temperature Tc 
occurs around the optimal doping, and then decreases in both the underdoped 
and overdoped regimes [T^; (b) a normal-state pseudogap metallic phase, where 

above Tc but 


14 
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an energy gap called the normal-state pseudogap Apg exists 
below the normal-state pseudogap crossover temperature T*. However, in contrast 


to the domelike shape of the doping dependence of Tc, T 
in the underdoped and optimally doped regimes [l4l-[^ 


is much larger than Tc 
and then monotonically 
decreases upon the increase of doping. In particular, measurements taken by us¬ 
ing a wide variety of techniques demonstrate that the normal-state pseudogap 
is present in both the spin and charge channels |14| ; (c) a normal-metal phase 
with largely transport properties. In the doped regime, charge carriers couple to 
spin excitations [^23 . The combined INS and resonant inelastic X-ray scattering 
(RIXS) experimental data have identified spin excitations with high intensity over 
a large part of moment space, and shown that spin excitations exist across the 
entire range of the SC dome [7,24 . However, the charge-carrier doping causes 
substantial changes to the low-energy spin excitation spectrum [^, while it has a 
more modest effect on the high-energy spin excitations 24 . In particular, RIXS 
experiments 24 show that the high-energy spin excitations persist well into the 


overdoped regime and bear a striking resemblance to those found in the parent 
compounds, indicating that a local-moment picture accounts for the observed spin 
excitations at elevated energies even up to the overdoped regime 24 . Experimen¬ 


tally, a large body of data available from a wide variety of measurement techniques 
have provided rather detailed information on cuprate superconductors, where some 
essential agreements have emerged. We refer the readers to the more detailed sum- 
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maries of experimental results available in the literatures 101 

Superconductivity, the dissipationless flow of electrical current, is a striking 
manifestation of a subtle form of quantum rigidity on the macroscopic scale j^, 
where a central question is what mechanism causes the loss of electrical resistance 
below Tc? It is commonly believed that the existence of electron Cooper pairs is 
the hallmark of superconductivity 31-33 , since these electron Cooper pairs be¬ 


have as effective bosons, and can form something analogous to a Bose condensate 
that flows without resistance. This follows from a fact that although electrons 
repel each other because of the Coulomb interaction, at low energies there can 
be an effective attraction that originates by the exchange of bosons [31] . In con¬ 
ventional superconductors, as explained by the Bardeen-Cooper-Schrieffer (BCS) 
theory 34,35 , these exchanged bosons are phonons that act like a bosonic glue 


to hold the electron pairs together, and then these electron pairs condense into a 
coherent macroscopic quantum state that is insensitive t o impu rities and imperfec¬ 
tions and hence conducts electricity without resistance 


34,351. The excitation in 


the SC-state has an energy gap A, which determines both the quasiparticle energy 
spectrum and the energy of the condensate |34H36l . In this conventional electron- 


phonon SC mechanism |35| , the resulting wave function for the pairs turns out to 
be peaked at zero separation of the electrons, which leads to that the SC-state has 
an s-wave symmetry [37| . As a consequence, the pairs in conventional superconduc¬ 
tors are always related to an increase in kinetic energy which is overcompensated 
by the lowering of potential energy 38 . At the temperature above Tc, the elec¬ 


tron is in the standard Landau Fermi-liquid state, which is generally referred to 
as a normal-state, where the de nsity of states near the Fermi level is smooth and 
generally treated as featureless 


34 


As in conventional superconductors, super¬ 
conductivity in cuprate superconductors results when charge carriers pair up into 
charge-carrier pairs, which is supported by many experimental evidences, includ¬ 
ing the factor of 2e occurring in the flux quantum and in the Josephson effect, 
as well as the electrodynamic and thermodynamic properties (M- 


45 . However, 


the normal-state of cuprate superconductors in the pseudogap phase is not normal 
at all, since the normal-state of cuprate superconductors in the pseudogap phase 
exhibits a number of the anomalous properties [14 -22 
not fit in with the standard Landau Fermi-liquid theory. 


in the sense that they do 
Superconductivity is an 


instability of the normal-state. However, one of the most striking dilemmas is that 
the SC coherence of quasiparticle peaks in cuprate superconductors is described 
by a standard BCS formalism, although the normal-state is undoubtedly not the 
standard Landau Fermi-liquid on which the conventional BCS electron-phonon SC 
mechanism is based. Angle-resolved photoemission spectroscopy (ARPES) experi¬ 
ments reveal sharp spectral peaks in the excitation spectrum |26[|27] , indicating the 
presence of quasiparticle-like states, which is also consistent with the long lifetime 
of electronic state as it has been determined by the conductivity measurements . 
Moreover, ARPES experiments also observe the Bogoliubov-type dispersion ofthe 
SC-state 46,^ predicted by the standard BCS formalism. However, as a natural 
consequence of the unconventional SC mechanism that is responsible for the high 
Tc, the ch arge-carrie r pairs in cuprate superconductors have a dominant d-wave 

This d-wave SC-state also implies that there is a strongly 


symmetry 39,42- 44 


Hi- 

3at- 


momentum-dependent attraction between charge carriers without phonons 
After over more than two decades of the painstaking effort, people are still de 
ing the very mechanism of superconductivity in cuprate superconductors, where 
the crucial issues in cuprate superconductors are (a) what is the nature of the glue 
binding charge carriers into charge-carrier pairs, such that they can travel macro¬ 
scopic distances without resistance? (b) whether the pseudogap has a competitive 
or collaborative role in engendering superconductivity? 

Very soon after the discovery of superconductivity in cuprate superconduc¬ 
tors 1^, Anderson [^ proposed that in the parent compounds of cuprate super¬ 
conductors, the spins form a superposition of singlets. This spin liquid of singlets 
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is so-called the resonating valence bond (RVB) state. Upon the charge-carrier 
doping, these RVB singlets would become charged, resulting in a SC-state. The 
RVB state is fundamentally different from the conventional Neel state in which the 
doped charge carrier can move freely among the RVB spin liquid and then a better 
compromise between the charge-carrier kinetic energy and spin exchange energy 


may be achieved. Since then many elaborations of this idea followed 48 -58 . In 


particular, it was realized that essential aspects of the RVB concept can be for¬ 
mulated within the charge-spin separation (CSS) slave-particle approach 5^62 

3 short 


since the essential physics of cuprate superconductors is dominated by the short- 
range repulsive interaction which remains relevant and causes CSS [6^. This 


CSS slave-particle approa ch also led to the development of the gauge theory for 
cuprate superconductors 


64-67 


However, in the framework of the original 
RVB theory 13, the AF e^^ange coupling J attracts electrons of opposite spins 
to be on neighboring sites. This is the result of states of very high energy with a 
spin gap, and the corresponding interaction has only high-energy dynamics 


The normal-state pseudogap is identified as the spin gap in the RVB state wrti 
an energy scale set by J, and therefore is associated with the breaking of the 
RVB singlets [^|^. In this case, Anderson suggested that the SC-state 
in cuprate superconductors is determined by the need to reduce the frustrated ki¬ 
netic energy of the system, where the strong frustration of the kinetic energy in the 
normal-state is partially relieved upon entering the SC-state, indicating that the 
kinetic energy causes superconductivity 


69 


On the other hand, it has been ar¬ 
gued phenomenologically that superconductivity in cuprate superconductors could 
arise from a lowering of the kinetic energy rather than the potential energy (70] . 
In this scenario, superconductors would exhibit qualitatively new features in their 
optical properties - a violation of the low-energy optical sum rule, and a change 
in the high-energy optical absorption when the system becomes SC 70,71 . Later, 


the high-precision optical measurements on cuprate superconductors in the near- 
infrared and visible region indicate small changes in the spectral weight associated 
with the onset of superconductivity |72||73 , therefore supporting this argument 70 
that changes in the kinetic energy are indeed occurring. In particular, the sim¬ 
ilar experimental results have by now been obtained by differently experimental 
groups |74|-l78|. More importantly, the recent experimental results IT^ from the 


ARPES measurements on cuprate superconductors indicate that Tc is correlated 
with the charge-carrier kinetic energy, which supports the notion of the kinetic 
energy driven superconductivity. Motivated by these experimental results 


72-78 


several calculations based on the strongly correlated models have been done to 
show that superconductivity may be driven by a lowering of the kinetic energy 
upon the formation of the SC-state 80-85 . By constructing an effective Hamil¬ 


tonian for spin polarons forming in weaMy doped antiferromagnets, it has been 
demonstrated that the driving mechanism which gives rise to superconductivity 
in such system is the reduction of the kinetic energy 80 . Moreover, a numer¬ 


ical study of the two-dimensional Hubbard model within the dynamical cluster 
approximation has shown the lowering of the kinetic energy below Tc for differ¬ 


ent doping levels 82 . These theoretical calculation 80 -85 shows that the paired 


charge carriers in the AFSRO background can be more mobile than the single 
charge carriers, and this can overcome the normal increase in the kinetic energy 
upon the pair formation 86 . 


Superconductivity in cuprate superconductors is something entirely new, a 
manifestation of the strong e lectron correlation, or Mottness [3^. In the early 


days of superconductivity, we [871 ®8] have developed a fermion-spin theory to con¬ 


front the strong electron correlation, where the constrained electron operator is 
decoupled as a product of a charge carrier and a localized spin, with the charge 
carrier represented the charge degree of freedom of the electron together with 
some effects of spin configuration rearrangements due to the presence of the doped 
charge carrier itself, while the spin operator represented the spin degree of freedom 
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of the electron, and then the strong electron correlation can be treated properly 
in actual calculations. In particular, these charge carriers and spin are gauge in¬ 
variant, and in this sense, the collective modes for these charge carriers and spins 
are real and can be interpreted as the physical excitations of the system. In this 
fermion-spin theory, the basic low-energy excitations are charge-carrier quasipar¬ 
ticles, the spin excitations, and the electron quasiparticles. In this case, the charge 
transport is mainly governed by the scattering of charge carriers due to spin fluc¬ 
tuations, and the scattering of spins due to charge-carrier fluctuations dominates 
the spin dynamics, while as a result of the charge-spin recombination, the electron 
quasiparticles are responsible for the electronic properties. Within the framework 
of the fermion-spin theory 87,88 , we have established a kinetic-energy driven 


SC mechanism 89-91 , where charge carriers are held together in d-wave pairs at 
low temperatures by the attractive interaction in the particle-particle channel that 
originates directly from the kinetic energy by the exchange of spin excitations in 
the higher powers of the doping concentration, and then these charge-carrier pairs 
(then electron Cooper pairs) condense to the d-wave SC-state. Although the phys¬ 
ical properties of cuprate superconductors in the normal-state are fundamentally 
different from these in the standard Landau Fermi-liquid state, the kinetic-energy 
driven SC-state still is conventional BCS-like with the d-wave symmetry, and then 
the obtained formalism for the charge-carrier pairing can be used to compute Tc 
and the related SC coherence of the low-energy excitations in cuprate supercon¬ 
ductors on the first-principles basis much as can be done for conventional super¬ 
conductors. Moreover, this kinetic-energy driven SC-state is controlled by both the 
charge-carrier pair gap and quasiparticle coherence, which leads to that Tc takes 
a domelike shape with the underdoped and overdoped regimes on each side of the 
optimal doping (^optimal ~ 0.15, where Tc reaches its maximum. On the other hand, 
the same charge-carrier interaction mediated by spin excitations that induces the 
SC-state in the particle-particle channel also generates the normal-state pseudo¬ 
gap state in the particle-hole channel [^. As a consequence, the SC gap and 
normal-state pseudogap coexist but compete in the whole SC dome. However, the 
normal-state pseudogap crossover temperature T* is much larger than Tc in the 
underdoped and optimally doped regimes, and then monotonically decreases upon 
the increase of doping, eventually disappearing together with superconductivity 
at the end of the SC dome. This kinetic-energy driven SC mechanism therefore 
provides a natural explanation of both the origin of the normal-state pseudogap 
state and pairing mechanism for superconductivity. 

It is beyond the scope of this article to provide an overview of various theo¬ 
ries of superconductivity in cuprate superconductors that have been put forward 
in the literatures, and some earlier reviews and different perspectives appear in 
Refs. (5 , [T5|[^ , 31,37,54,58,62,67,93,94 . In this article, we only attempt to review 
compr^ensively the general framework of the kinetic-energy driven SC mechanism 
in the context of our work [87| - [92| and to summarize several calculated results of 
physical quantities obtained based on the kinetic-energy driven SC mechanism. 
The number of topics for this review article is listed as follows. In section 
we |87l[88 give an overview of the fermion-spin theory, and show that within the 


decoupling scheme, the fermion-spin representation is a natural representation of 
the constrained electron defined in a restricted Hilbert space withou t double elec¬ 
ts introduced 


tron occupancy. The kinetic-energy driven SC mechanism [89 -91 
in section In the rest of sections, we show how this kinetic-energy driven SC 


mechanism yields many results that are in broad agreement with various key ex¬ 
perimental facts observed on cuprate superconductors. Superconductors are not 
only perfect conductors, but also exhibit the so-called Meissner effect, where they 
expel magnetic fields. In section we consider main features of the doping de¬ 
pendence of the electromagnetic response observed on cuprate superconductors by 
using the muon-spin-rotation measurement technique, and show that in analogy to 
the domelike shape of the doping dependence of Tc, the maximal superfluid density 
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Ps occurs around the critical doping (^critical ~ 0.195, and then decreases in both 
lower doped and higher doped regimes. In section we turn to the comparison of 
the calculated result of the dynamical spin response with RIXS-INS experimental 
data. It is shown that the low-energy spin excitations in the SC-state have an 
hour-glass-shaped dispersion, with commensurate resonance that appears in the 
SC-state only, while the low-energy incommensurate (IC) spin fluctuations can 
persist into the normal-state. The high-energy spin excitations in the SC-state 
on the other hand retain roughly constant energy as a function of doping, with 
spectral weights and dispersion relations comparable to those found in the par¬ 
ent compounds. A brief description of the interplay between superconductivity 
and normal-state pseudogap state is given in section where we identify the 
normal-state pseudogap as being a region of the self-energy effect in the particle- 
hole channel in which the normal-state pseudogap suppresses the spectral weight 
of the low-energy excitation spectrum. This normal-state pseudogap disappears at 
T*, and then system crossovers to the normal-metal phase with largely transport 
properties at the temperatures T > T*. In section we discuss the effect of the 
normal-state pseudogap on the infrared response of cuprate superconductors, and 
show that in the underdoped and optimally doped regimes, the transfer of the part 
of the low-energy spectral weight of the conductivity spectrum to the higher energy 
region to form a midinfrared band is intrinsically associated with the emergence 
of the normal-state pseudogap. Finally, the article concludes with the suggestions, 
in section for future work. 

2. Fermion-spin theory 
2.1. Model 

In cuprate superconductors, the single common feature in the layered crystal struc¬ 
ture is the presence of one to several Cu 02 planes in the unit cell [iH^, and it 
seems evident that the nonconventional behaviors of cuprate superconductors are 
dominated by the Cu02 plane. In this case, as originally emphasized by Ander¬ 
son [3 , the essential physics of the doped Cu 02 plane is contained in the one-band 
larged Hubbard model on a square lattice, 

H = - ^2 M X] ^ X] (1) 

<IV>(T Icr I 

where the summation is over all sites I, and the hopping integrals connect sites 
I and I'. We will restrict to our attention to the nearest (t) and next nearest {—t') 
neighbor hopping. Cj^ and Cia are electron operators that respectively create and 

annihilate electrons with spin a, nia = Cj^Ci^, and p is the chemical potential. 
This large-C/ Hubbard model Q indicates that the interactions in cuprate super¬ 
conductors are dominated by me on-site Mott-Hubbard term U, which is very 
large as compared with the electron hopping integrals t and t', i.e., U ^ t, t', 
and therefore leads to that electrons become strongly correlated to avoid double 
occupancy. In this case, the on-site Mott-Hubbard term must be dealt properly 
before bothering with relatively minor terms (6^|96j . It has been shown that 
the correct way to deal with this large-C/ term in Eq. 0 is to renormalize it by 
means of a canonical transformation e®'^, which eliminates large-C/ term from the 
block which contains no doubly occupied states, and which presumably contains 
all the low-energy eigenstates and thus the ground state, and then the transformed 
Hamiltonian can be obtained as, 

^ E ClCi+^. + t'Y, ClCi+,^ + pY1 E S' • (2) 

IfjG Ira la Iff 
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with the nearest-neighbors fj = ±x, ±y, the next nearest-neighbors f = ±x zb y, 
the magnetic exchange coupling constant J = 41“^jU, the spin operators S; = 
{Sf,Sj,Sf). The kinetic-energy term in Eq. Q describes mobile charge carriers 
in the AF background, while the Heisenberg term in Eq. Q describes AF coupling 
between localized spins. In particular, the nearest-nei^bor hopping integral t 
in the kinetic-energy term is much larger than the magnetic exchange coupling 
constant J in the Heisenberg term, and therefore the spin configuration is strongly 
rearranged due to the effect of the charge-carrier hopping t on the spins, which 
leads to strong coupling between the charge and spin degrees of freedom of the 
electron. This transformed Hamiltonian (§ is so-called t-J model acting on a 
restricted Hilbert space without double electron occupancy, where there are three 
physical states only, 


|0 >, I t>, I i> . (3) 

The essential physics is cuprate superconductors obeying the t-J model Q, in 
which the hopping integrals of t and t' to states having higher energy U are removed 
in favour of the magnetic exchange interaction J [96] . At half-filling, this t-J model 
Q is reduced to an AF Heisenberg model, where the degree of freedom is local 
spin only. In particular, it has been demonstrated explicitly the local SU(2) gauge 
invariance of the Heisenberg model written in terms of electron operators with a 
constraint of one particle per site 1^. In spite of its simple form, the t-J model Q 
has been proved to be very difficuTFto analyze, analytically as well as numerically, 
because of the restriction of the motion of electrons in the restricted Hilbert space 
without double electron occupancy. 


2.2. Constraints and sum rules 


The strong electron correlation originates from a large on-site repulsion between 
two electrons occupying the same site in the Hubbard model Q, which effectively 
translates into an elimination of double occupancy in the t-J model ([^. There 
are two ways to implement the crucial requirement of no double occupancy (87| : 
either to solve the t-J model Q combined with a single occupancy local constraint 


60-62 


E 


clcu, < 1 , 


(4) 


or to introduce the constrained electron operator replacing C;o- by, 

Ci„ = Ci„{l-ni_„). (5) 

In this section, we will use both representations to clarify this matter. As a 
consequence of the electron motion in the restricted Hilbert space ^ without 
double occupancy, the constrained electron operator satisfies following relation 
|87l, 




( 6 ) 


this leads to a sum rule for the constrained electron, 

(E Alo.) = 1-4. 


(7) 


where 5 is the charge carrier doping concentration, and (• • •) means thermody¬ 
namical average. On the other hand, the constrained electron operator Cia obeys 
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a special on-site anticommutation relation, 


= (8) 

CT <J 

with its expectation value, 




cl}) = 1 + 5, 


(9) 


which gives rise to a sum rule for the electron spectral function ^o-(k, cu) 


E 


du 


Ao-(k, Lu) — 1 -|- (5. 


( 10 ) 


In the Hubbard model Q, the large-t/ forces a finite density of the many-electron 
states to split out of the band continuum of states to the high-energy side, forming 
the upper Hubbard band 96 , then the remaining continuum (lower Hubbard band) 
is completely described by the t-J model Q. In this case, a doped charge carrier 
leaves behind an empty site. However, each such empty site can be occupied only 
by either a spin-up or spin-down electron in the restricted Hilbert space Q. In 
particular, just from the charge-carrier doping, the empty part of the electron 
spectrum at low energies has a weight of 26. Furthermore, as the charge-carrier 
doping annihilates one state in the filled part of the electron spectrum, there are 
1 — 6 electron state (per site) remaining below the Fermi surface. This gives rise 
to a total weight of the lower Hubbard band of 1 -|- 5 [^[^. Eqs. 0 and @ 
are the exact sum rules for the t-J model 0, and they ^ould be preserved in an 
adequate treatment. 


2.3. Slave-particle theory 

The high complexity in the t-J model 0 comes mainly from the electron single 
occupancy local constraint S- This electron single occupancy local constraint 
0 can be exactly taken into acco unt only by n umerical methods, such as the 


variational Monte Carlo technique 48,49,51,58, exact cluster diagonalization 


(Too , and various realizations of the quantum Monte Carlo method jlOl]. However, 


the exact diagonalization is limited by system sizes, while the quantum Monte 
Carlo technique faces the negative sign problem for lower temperatures. Apart 
from these numerical techniques, an intuitively appealing approach to implement 
this electron single occupancy local constraint 0 and the CSS scheme is the slave- 
particle approach (^-62 , where the physics ofno double occupancy is taken into 


account by representing the constrained electron as a composite object created by. 


Cia = a]fi, 


I J la •, 


( 11 ) 


aj as the fermion 


with aj as the slave boson and fia as the fermion or vice versa, i.e., 
and fia as the boson. This way the nonholonomic constraint 0 is converted into 
a holonomic one 62 


+ = ( 12 ) 

a 

which means a given site cannot be occupied by more than one particle. In this 
slave-particle representation the charge degree of freedom of the constrained 
electron is described by the operator ai, while the spin degree of freedom of the 
constrained electron is described by the operator //o-, and then the elementary 
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charge and spin excitations are so-called holon and spinon, respectively. However, 
a new ?7(1) gauge degree of freedom must be introduced to incorporate the single 
occupancy local constraint (12), which means that the slave-particle representation 
should be invariant under a local U{1) gauge transformation, 


ai 




idl 


(13) 


and then all physical quantities should be invariant with respect to this transforma¬ 
tion. This reflects that the holon a; or spinon itself is not gauge invariant, and 
they are strongly coupled by the 17(1) gauge field fluctuation [^,^65 67 


In this 
therefore they 


Moreover, 


sense, the collective modes for the holon and spinon are not real arr_ 

can not be interpreted as the physical excitations of the system [l02 
there are a number of difficulties in this slave-particle approach. First of a ll, in the 
slave-boson version, the AFLRO correlation is absent for zero doping [103 


so that 


the ground-state energy is high compared with the numerical estimate 49,104 
and the Marshall sign rule 105] is not obeyed. Alternatively, in the slave-fermion 


approach, the ground-state is the stat e with AFLRO for the undoped case and per¬ 
sists until very high dopin g ('^ 60 %) |106|. In particular, it should be noted that 


in the actual calculations 61,62], the electron single occupancy local constraint 


(12) is explicitly replaced by a^obal constraint, and therefore the representation 
space is much larger than the restricted Hilbert space ([^ for the physical electron. 
The local nature of the constra int is of p rime importance, and its violation may 
lead to some unphysical results [107 -109 . This is why the crucial requirement for 
the t-J model is to impose the electron single occupancy local constraint (|^. 


2.4. CP^ representation 

For convenience in the following discussions at this section, another useful approach 
to implement the elec tron single occupancy local constraint ([^ and the CSS scheme 


has been developed [no 
decoupled as 


where the physical electron in the t-J model ([^ 


IS 


Cia = h\h^, (14) 

supplemented by the local constraint, 

= (15) 

cr 

with the spinless fermion hi keeping track of the charge degree of freedom of the 
constrained electron (holon), while the spinful boson keeping track of the spin 
degree of freedom of the constrained electron (spinon), and then the electron single 
occupancy local constraint in Eq. (|^, 

Y,ClCi^ = l-h]hi<l, (16) 


is satisfied, where = h\hi is the holon number at site I, equal to 1 or 0. This 
decoupli ng s cheme, the so-called slave-fermion CP^ representation, was proposed 
in Ref. 110 . However, it is similar to the slave-particle representation ( jll] ), the 


spinons and holons in the slave-fermion CP^ representation are not physical ob¬ 
jects, and there is an arbitrary phas e rel ation between the two related to the single 
occupancy local constraint (15) as [ITo] , 


hi hie^^h bia —?• biaC 


idl 


(17) 
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This local U{1) gauge degree of freedom, of course, is canceled for the physical 
electron. However, it should be noted [STIISS] from Eq. (15) that so long as 


h\h, = hlh,, = l,Y.,clQ, = Q. no matter what is the valne In the 

slave-fermion CP^ representation, the choice Yla = 1 is convenient, because 
it also guarantees the local condition = 1; when h\hi = 0. However, 

the local constraint = 1 means the presence of one boson (spin-up or 

down) on each site, i.e., a spin even to an empty site has been assigned. This 
will not affect the physical expectation values, because the charge carrier number 
expectation {h\hi) = {h\^hi„) will remove all spurious effects. Never thel ess, the 


extra degrees of freedom will affect the partition function. In subsection 2.7 we 
will define a projection operator to cure this defect. As a result, the sum rules 
and (10) will be satisfied exactly. 


2.5. Fermion-spin transformation implement gauge invariant charge carrier 

The decoupling of the charge and spin degrees of freedom of ele ctron is undoubt¬ 
edly correct in one-dimensional interacting electron systems 111, where the charge 
and spin degrees of freedom of electron are represented by boson operators that 
describe the excitations of charge-density wave and spin-density wave, respectively. 
In particular, the typical behavior of the non-Fermi-liquid, showing up as CSS and 
vanishing of the quasiparticle re sidu e, has been demonstrated theoretically within 
the one-dimensional t-J model 112|. Moreover, the excitations of the charge and 


spin degrees of freedom of electron as the real elementary excitations in the one¬ 
dimensional cuprates has been observed directly by the ARPES experiment [113 


Therefore both theoretical and experimental studies indicate that the existence of 
the real excitations of the charge and spin degrees o f freedom of electron is com¬ 
mon in one-dimensional interacting electron systems [111 -114 . However, the case 
in cuprate superconductors (two-dimensional strongly correlated electron systems) 
is very complex. Among the anomalous properties of cuprate superconductors in 
the normal-state in the underdoped and optimally doped regimes, a hallmark is 
the charge transport [^[6,14 , where the low-energy conductivity deviates strongly 


from the Drude behavior, and is carried by 6 charge carriers. This is very natural 
in that the low-energy spectral weight of the conductivity must vanish when 5 —)• 0. 
It follows that a superconductor that forms out of the under doped cu prates must 
have a superffuid density ps given by this spectral weight 13,115,116 , so that the 
superfluid density ps in the underdoped regime vanishes more or less linearly with 
the decrease of the charge-carrier doping concentration 6. This in turn gives rise 
to the linear relation b etwee n and ps observed in cuprate superconductors in 
the underdoped regime 116 . In corresponding to the non-Drude behavior of the 


low-energy conductivity in the underdoped and optimally doped regimes, the re¬ 
sistivity exhibits a linear temperature behavior over a wide range of temperatures. 
It has been argued that these experimental facts are a strong experimental evi¬ 
dence snpporting the notion of CSS, since not even conventional electron-electron 
scattering would show the striking linear rise of scattering rate above the Debye 
frequency, and if there is no CSS, the phonons should affect these properties [63] . 
In this case, a formal theory with the gauge invariant excitations of the charge 
and spin degrees of freedom of electron in the two-dimensional strongly correlated 
electron systems, i.e., the issue of whether the excitations of the charge and spin 
degrees o f freedo m of electron are real, is centrally important |l02| . In this subsec 
start from the above slave-fermion CP^ 


tion, we 87,88 


approac. 


(|14|) and show 
y, the con- 


that if the electron single occupancy local constraint is treated proper 
strained electron in the t-J model ([^ can be decoupled by introducing the charge 
carrier and spin, where the collective mode for the charge carrier or spin is real 
and therefore can be interpreted as the physical excitation of the system. 

First of all, we |87j examine the properties of the spinless bosons 6/. The spinless 
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boson creation and annihilation operators are expressed in the infinite-dimensional 
Fock space as 117|, 


b] = 


( 0 
^/T 
0 
0 

V - 


0 

0 

0 


0 

0 

0 

V3 


\ 


/ 


hi = 


/ 0 \/T 0 0 ... \ 

0 0 \/2 0 

0 0 0 \/3 

0 0 0 0 

V./ 


(18) 


However, if the boson occupation number is restricted to 0 or 1, the spinless boson 
creation and annihilation operators in Eq. (18) are reduced immediately to the 


two-dimensional space as 87 


^\ = 


0 0 
1 0 


hi = 


(19) 


which are nothing but the spin-lowering and spin-raising operators for spin 
S = 1 / 2 . These spin operators and behave as fermions on the same site, 

and as bosons on different sites, and therefore satisfy the hard-core constraints 
hbi = b\h\ = 0. 

Now we turn to explore further the properties of the CP^ bosons 6 ;^- Since 
the CP^ bosons bi^ and 64 satisfy the CP^ local constraint Yla = 1 in Eq. 

(15), the empty and doubly occupied states have been ruled out, and only the 


spin-up and spin-down singly occupied spin states are allowed. In particular, due 
to the symmetry of the spin-up and spin-down states, | occupied)'|' = f q 


and 


em 

0 

1 


pty)t ~ 1 ^ singly-occupied and empty spin-up, while | occupied)^^ = 


and I empty) 4 . = 


are singly-occupied and empty spin-down states, 


local constraint (15) can be represented in the basis 




bi = I occupied)'!' 4 .(occupied |= e®'*’-'- 


occupied)^^ -j-(occupied |= e®®"’" ^ 


2 $ I 


1 
0 

0 0 
1 0 

0 1 
0 0 


and 




and then all the hard-core boson conditions, i.e. 
biaha = 0, (without summation over a), are satisfied. As a result, the CP^ boson 




+ b\j)la 



87 

88 

75 -, 

( 20 a) 


( 20 b) 

= 1 , b] b] 

’ 1(7 2 

tJ 


operators 64 and 64 together with the local constraint (15) are identified with 


the spin lowering operator with an additional phase factor e®'*’*'!' and raising 
operator with an additional phase factor e®^'-i-, respectively. Consequently, the 
corresponding CP^ ordinary boson occupation space together with the local con¬ 
straint of one boson per site is identified with the natural spin 1 /2 representation 
space, while the phase factor in Eq. (20) is closely related to this transfor¬ 
mation of the representation spaces, and therefore carries some messages of the 
spin degree of freedom of the constrained electron, especially, some messages of 
the spin configuration rearrangements in the doped case [^. In this case, the 
electron decoupling form in Eq. (14) can be expressed as. 


C4 = h\e^^^^Sr 


i > 


Cii = h\e^^ 


4S+, 


( 21 ) 
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while the local ?7(1) gauge transformation © therefore is rewritten as, 

hi —T” hie^^‘, —)• + 9i. 


( 22 ) 


In particular, the phase factor in Eq. (21) can be incorporated into the 

spinless fermion operator h], and then as the solution of the electron single oc¬ 
cupancy constraint g E.ClCi^ < 1 under CPi sla^^ermion convention ^ 


So- = 1, we find the following transformation 


87 


88 


Of = dtsr. c,i = his*, 


(23) 


where the spinful fermion operator hi^ = hi represents the charge degree of 

freedom of the constrained electron together with some messages of the spin degree 
of freedom (charge carrier), while the spin operator Si represents the spin degree of 
freedom of the constrained electron, and then the electron single occupancy local 
constraint 0 . 

= s+'«TdtSr+sr'>4'>lE = '>i'4(stsr+srs,+) 

<J 

= l-h\hi<l, (24) 

is always satisfied in actual calculations. In other words, the electron single 
occupancy local constraint (|^ is implemented exactly using the transformation 
(23). _TI^ electron decoupling form (23) is called the fermion-spin transforma- 

obey the 


tion 


87 


88 


Since the spinless fermion hi and spin operators and Sf 

anticomiinutation relation and Pauli spin algebra, respectively, it is then easy to 
show that the spinful fermion hi^j also obeys the same anticommutation relation 
as the spinless fermion hj. However, in contrast to the holon and spinon in the 


slave-particle approach (11) and the slave-fermion CP^ formalism (14), the charge 
carrier hia or spin Si itself is invariant under the local U{1) gauge transformation 
m, i.e.. 


hia = hie 




hie^^^ 




= hie-^'^^^ = hi„, Si ^ Si, 


(25) 


which leads to that all physical quantities from charge carriers or spins are invariant 
with respect to the gauge transformation (22). In this sense, the collective mode 
for the charge carrier or s pin i s real and therefore can be interpreted as the physical 
excitation of the system 


102 


The essential physics of the fermion-spin transformation ( |23[ ) is simple: at half¬ 
filling, the t-J model (1^ is rednced as an AF Heisenberg model) where each lattice 
site is singly occnpied by a spin-up or spin-down electron, then the electron spins 
(the spin degree of freedom) are coupled antiferromagnetically with AFLRO. With 
the charge-carrier doping, the spin configuration must be rearranged to provide 
the lowest magnetic energy. However, at the same time, the holon (the charge 
degree of freedom) can feel self-consistently the effect of the spin configuration 
rearrangements through the strong coupling between the two degrees of freedom, 
and then the charge carrier arises from the binding of the holon and the phase 
factor carrying some messages of the spin configuration rearrangements due 
to the presence of the doped charge carrier itself. In particular, when a charge 
carrier at site I hops to site I' in the t-J model ([^, the spin configurations are 
optimized by the fermion-spin transformation ( |23[ ) , indicating that the frustration 
effect of charge carrier hopping in the spin background is well taken care of by 
the fermion-spin transformation ( pS] ). In fact, the representation of the hard-core 
boson in terms of spin raising and lowering operators is essential, because whenever 
a charge-carrier hops it gives rise immediately to a change of the spin background 
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as a result of careful treatment of the electron single occupancy l ocal c onstraint. 
This is why the t-term is so efficient in destroying the AFLRO [^, 118 


The fermion-spin theory (23) also indicates that the constrained electron is a 


composite object, with charge carriers and spins being the physical excitations. 
However, the charge-carrier quasiparticle and spin excitation are strongly renor¬ 
malized each other because of the coupling between the two degrees of freedom. In 
this case, three basic low-energy excitations for the charge-carrier quasiparticles, 
the spin excitations, and the electron quasiparticles, respectively, emerge as the 
propagating modes in a doped Mott insulator [^, with the charge-carrier quasi¬ 
particles that are responsible for the charge transport, and the spin excitations 
dominate the spin response, while as a result of the char ge-spin re combination, 
the electron quasiparticles govern the electronic properties (TTqHT^ . 


2.6. Fermion-spin representation - a natural representation for constrained elec¬ 
tron 


Now we show that in the decoupling scheme^ the fermion-spin transformation (23) 


is a natural representation for the constrained elec tron defined in a restricted 
Hilbert space without double electron occupancy 122 . In Eq. © , the constrained 
electron operators and Cia are expressed in terms of the unconstrained electron 
operators and Cia as (5/^ = — and Qo- = — respectively. 

Although the constrained electron operator {Cia) does not create (destroy) any 
doubly occupied sites 96 , the unconstrained electron operators and C/o- are 


thought to be operators to operating within the full Hilbert space. In particular, 
the constrained electron operators Cq and C 4 can be rewritten as 


96 


Cii = Ci^{l-ni^) = Ci^C\^Ci^ = Ci^Sl, 


(26a) 

(26b) 


where the spin index | (f) of the unconstrained electron operator Cq (Cq) in 
the right-hand side in Eq. (26) is not an independent degree of freedom, since 
the spin fluctuation of the system is mainly described by the spin operator Si. In 


other words, in the constrained electron operators (26 
operator Cia in the right-hand side in Eq. (26) main 


the unconstrained electron 
y describes charge degree of 
freedom of the constrained electron together with some messages of the spin degree 
of freedom, while the spin operator Si represents the spin degree of freedom of the 
constrained electron. Furthermore, the constrained electron operators in Eq. (26) 
are exactly same as quoted in Eq. (23) in the fermion-spin representation if the 
constrained electron (26) is decoupledaccording to its charge and spin degrees of 
freedom, i.e., in the decoupling scheme, the unconstrained operator Cia and spin 
operator Si in the right-hand side in Eq. (26) are decoupled as two independent 
operators Cia and Si , and then Cia and Si in Eq. (26) commute each other in the 
decoupling scheme. To see this point clearly, the constrained electron operators 


Cq and Cq in Eq. (26) in the decoupling scheme can be rewritten in terms of a 
particle-hole transformation for the unconstrained operator Ci„ —?• h\_^ as. 


Of = dtsr. c,i = fct s+ 


(27) 


which are exactly same as quoted in Eq. (23) in the CSS fermion-spin representa¬ 
tion. 


2.7. Projection operator 

In the local representation, the restricted Hilbert space without double electron 
occupancy in the t-J model ([^ consists of three states, |0), | f)) I i) as we have 
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mentioned in Eq. Q. However, in the fermion-spin transformation (23), there 
are four states [charge) 0 [spin), namely |l,t)) |0)t)) and |0,4,), wEere 1 or 

0 means charge-carrier occupation or empty. In this case, a pr oje ction operator 
P can be introduced to remove the extra degrees of freedom 1^. The matrix 
elements of this projection operator can be defined as. 


Pna ^ \l<-){a\, (28) 

where \k) is one of the bases of the physical states, while |a) is one of the bases in 
the fermion-spin representation space [charge)® [spin). Since the space dimensions 
of [k) and [a) are different, the usual relations for the projection operator P^ = 
P = pi are not satisfied. Using this projection operator, the constrained electron 
operators in the restricted Hilbert space of no double occupancy can be defined 
as, 


C„ = FifttSfP/. Cl^ = P,h„S+p!, 
Qi = p,h],s+p!, c* =p,ft4srpt. 


(29a) 

(29b) 


where Pi is the projection operator for the site I, and P^ is the Hermitian conjugate 
of Pi- Making use of the matrix representation of the charge-carrier operators. 


hl = 


0 1 
0 0 




0 0 
1 0 


(30) 


all these oper ators (29) in matrix form can be obtained explicitly (see the Appendix 
Appendix A). In particular, the constrained electron operators in the basis of the 
physical states [0), [ f), [ i) can be written as [87], 



pi — 


pi — 

'-"4 “ 



(31a) 

(31b) 


which are nothing but the Hubbard X operators ^ 04 .) dc. |95| . It is then 
straightforward to check that 


T.6IP 


= 1 — n 


(h) 

I ’ 


E(Cl».C,l} = l + n 


{h) 

I ’ 


(32a) 

(32b) 


with the charge-carrier number operator, 


1 0 0 
=[000 
0 0 0 


(33) 


Taking expectation values of Eq. (32), one sees immediately that the sum rules 0 
and (10) are exactly satisfied. It is therefore shown that the fermion-spin trans¬ 


formation (23) defined with an additional projection operator P satisfies exactly 


the single occupancy local constraint and all sum rules, i.e., they are an exact 
mapping (^. 
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However, the projection operator P is cumbersome to handle in many cases, 
and it has been dropped in the actual calculations (OT-92 . Now let us see which 
of these properties are still preserved and what kind of errors we are committing 
in such approximate treatments. First of all, the electron single occupancy local 
constraint (24) is exactly obeyed even in the mean-field (MF) level 87, [^. For 
comparison, it should be noted that in the usual slave-particle approac 


103 |106 , the electron single occupancy local constraint is explicitly replacedby a 



global constraint in the MF approximation, and therefore the representation space 
is much larger than the representation space for the physical electron. From this 
point of view, the treatment of the electron single occupancy local constraint for the 
physical electron in the fermion-spin theory is m uch b etter than the slave-particle 
approach, and therefore the MF results 87,88 120 based on the fermion-spin 
theory even without projection operator are better than those obtained within the 
framework of the slave-particle MF theory. Secondly, those expectation values of 
electron operators, including spin-spin correlation functions, which should vanish, 
actually do not appear due to the presence of the charge-carrier number operator 
= hj^hio- = hjhi- By adding the extra spin degree of freedom to an empty 
site, we are making errors of the order 5 in counting the number of the spin states 
(8711^1 [ 12 ^, which is negligible for small doping. To show this point clearly, we 
can map electron operator Qo- with the electron single occupancy local constraint 

supplemented by the 


onto the slave-fermion formulism 


62 


as Ci(j aiQ 


local constraint h]hi + ~ solve the local constraint in the 

slave-fermion formulism by rewriting the boson operator a/o- in terms of the CP^ 

boson operator as ai„ = — h\hi supplemented by the local constraint 

'^^b\j)i„ = 1 . As we have mentioned in subsection 2.5, the CP^ boson operators 


hi^ and 64 with the local constraint = 1 can be identified with the 

spin lowering and raising operators, respectively, defined with an additional phase 
factor, therefore the projection operator is approximately related to the charge- 

carrier number operator by Pi ~ — h\^hia = 1 — h\hi, and its main role is 

to remove the spurious spin when there is a charge carrier at a given site I [ 88 ] . 
Thirdly, the essential physics of the local U{1) gauge invariant charge carrier and 
spin is still kept 88 , which guarantees that the collective mode for the charge 


carrier or spin is r eal and therefore can be interpreted as the physical excitation 
of the system [102 . 


2.8. Summary 

In this section, a brief review is given for the fermion-spin theory. In the framework 
of the fermion-spin theory, the constrained electron is decoupled as a product of a 
gauge invariant charge carrier and a localized spin, with the charge carrier keeping 
track of the charge degree of freedom together with some effects of the spin config¬ 
uration rearrangements due to the presence of the doped charge carrier itself, while 
the spin operator keeping track of the spin degree of freedom. The main advantage 
of this fermion-spin theory is that the electron on-site local constraint for single 
occupancy is exactly satisfied even in the MF approximation. In particular, in 
the decoupling scheme, the fermion-spin representation is a natural representation 
of the constrained electron defined in a restricted Hilbert space without double 
electron occupancy. 


3. Kinetic-energy driven supercondncting mechanism 

The theory of superconductivity in cuprate superconductors remains one of the 
most important outstanding problem in the field of condensed matter physics. 
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From the experimental front, the measured result from optical spectroscopy on 
cuprate superconductors with simultaneous time and frequency resolution demon¬ 
strates that bosonic excitations o f ele ctronic origin are the most important factor 
in the formation of the SC-state 124| . In particular, the combined RIXS-INS ex¬ 
perimental data have identified spin excitations with high intensity over a large 
part of momentum space, and shown that spin excitations exist across the en¬ 
tire range of the SC dome (7,24 , which provides a clear link between the pairing 
mechanism and spin excitations. In this case, a question is raised whether the spin 
excitation, which is a generic consequence of the strong Coulomb interaction in the 
large-C/ Hubbard model Q, can mediate charge-carrier pairing in cuprate super¬ 
conductors in analo gy to me phonon-mediate pairing mechanism in conventional 
superconductors 31 125 ? On the theoretical side, it has been shown that the BCS 
theory is not specific to a phonon-medi ated int eraction, and other bosonic excita¬ 
tion can also serve as the pairing glue 


31,125 


In this section, we review briefly 
the kinetic-energy driven superconductivity |89| - [^ , where the effective attractive 
interaction between charge carriers originates in their coupling to spin excitations, 
and therefore the spin excitation has been served as the pairing glue. 

Since the t-J model ® is obtained from the Hubbard model Q as mentioned 
in section taking the large-t/ limit and making certain approximations (^ , 
there is a mixing of kinetic energy and potential energWCoulombic contribution) 
in going from the Hubbard model Q to the t-J model Q, i.e., the original kinetic 
energy in the Hubbard model Q nas been reorganizea as the kinetic energy of 
the lower Hubbard band in the t-J model which therefore contains a strong 
Coulombic contribution due to the restriction of no doubly occupancy of a given 
site. On the other hand, the magnetic exchange energy in the t-J model Q also 
has mixed electronic kinetic and Coulombic origins, since the intersite magnetic 
exchange interaction in the t-J model Q is due to the mixing of states in the 
lower and upper Hubbard bands on adjacent sites, i.e., the virtual hopping, which 
is enabled by AF alignment, and involves kinetic energy |86|. In the fermion-spin 
representation (fls]), the t-J model Q can be rewritten as 




l-\-r 


Ifj 


It 


+ A;+,,A4srs+,) - fE '■L'*'- + •'=11E ■ s. 


I+VI 


(34) 


li) 


where Jeff = (1 — 5)^J, 6 = {hj^hi„) = {h]hi) is the charge-carrier doping con¬ 
centration, while the electron degrees of freedom in the original t-J model (1^ has 
been replaced by two degrees of freedom associated with the doped charge car¬ 
rier within a spin correlation background related to the electronic structure of the 
AF correlation phase. As an important consequence, the mixing of kinetic energy 
and potential energy in the kinetic-energy term of the original t-J model ([^ has 
been released as the interaction between charge carriers and spins in the t-J model 
(34), which therefore dominates the essential physics of cuprate superconductors, 
while the magnetic exchange energy term is to form an adequate spin configuration 
only [63[|69] , which is strongly rearranged due to the effect of the charge-carrier 
hopping t on the spins. In particular, since the electron single occupancy local con¬ 
straint (then the strong electron correlation) in the t-J model (34) has been taken 
into account properly within the fermion-spin representation (23), the remaining 
analysis, including the treatment of the interaction between charge carriers and 
spins in the t-J model (34), can be fully carried out by a perturbation theory 
based on the conventional many-body technique as we do for the electron-phonon 
interaction underlying conventional superconductors. 

In a superconductor, electrons form Cooper pairs with an energy gap in the 
single-particle excitation spectrum when they become superconductors. This gap 


U 
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corresponds to the energy for breaking a Cooper pair of electrons and creating two 
quasiparticles, therefore measures the strength of the binding of two electrons into 
a Cooper pair. In particular, electron Cooper pairs in conventional superconduc¬ 
tors are formed in momentum space, not in real space, because the onset of the 
long-range phase coh erenc e occurs due to the overlap of the electron Cooper pair 
wavefunctions [M, 35 126 . As a consequence, the order parameter in the charac¬ 


terization of the SCtransition and the Cooper pair wavefunctions are the same: 


the order parameter is a magnified version of the Cooper pair wavefunctio ns |126 
On the other hand, within the kinetic-energy driven SC mechanism |89f|91] , the 
charge-carrier pairing also occurs in momentum space, then the order parameter 
in the characterization of the SC transition in cuprate superconductors and the 
charge-carrier pair wavefunctions are the same, i.e., the charge-carrier pairing and 
the onset of the phase coherence take place simultaneously at the charge-carrier 
pair transition temperature, then the SC transition temperature Tc is identical to 
the charge-carrier pair transition temperature. However, cuprate superconductors 
are doped Mott insulators, where charge carriers (then electrons) couple to spin 
excitations [7 


23 . In this case, we first consider electron Cooper pair in real space. 


since it can give some insights into the nature of the SC-state. For example, the 
amplitude of the electron Cooper pair between the elec trons a t sites I and I' can 
be expressed in the fermion-spin representation (23) as 


89-91 


A(( - 0 = - CjC/,,) = {h,-,h,,iS+S^ - fc4VtSrS,t). 


''t r*! \ _ 


(35) 


In the doped regime without an AFLRO, charge carriers move in the spin liquid 
background, where the spin correlation functions (Sfi 5 ,7) = = xi-i': and 

then the amplitude of the electron Cooper pair in Eq. (|35|) can be rewritten as, 


A(/-0 = -M-z'Ah(/-0, 


(36) 


with the amplitude of the charge-carrier pair, 

Ah(^ — I') = — hii^hifi), 


(37) 


which shows clearly that the amplitude of the electron Cooper pair is closely 
related to the corresponding one of the charge-carrier pair, and is proportional to 
the number of the doped charge carriers, but not to the number of electrons. In this 
case, the essential physics in the SC-state of cuprate superconductors is dominated 
by the corresponding one in the charge-carrier pairing state. The charge-carrier 
pairing in momentum space can be c onsidered as a collective phenomenon, while 
that in real space as individual [126 . Of course, the transition into the SC-state 
always occurs in momentum space. However, in the extremely low-doped regime 
with an AFLRO, where the spin correlation functions Sfi) {Sf S'jt), and then 
the conduct is disrupted by this AFLRO. In the following discussions, we 
only focus on the case without AFLRO, and leave the case with AFLRO for 
further discussions. 


91 

the 


3.1. Charge-carrier and spin Green’s functions 


For convenience in the following discussions, we first define the charge-carrier di¬ 
agonal and off-diagonal Green’s functions as [88], 


g{l -l',t- t') 

r(/-/',t-t') 

T\l-l',t-t') 


{{hia{t);hl^{t'))), 
{{hiiit); hi>^{t'))), 

((/.;,(*); (('))), 


(38a) 

(38b) 

(38c) 
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respectively, and the spin Green’s functions as, 

= {{S+{ty,Si;it'))), (39a) 

= {{Sf{t);S!,{t'))), (39b) 


where ((.. .)) is an average over the ensemble. Since the spin system in the t-J 
model (34) is anisotropic away from half-filling, therefore two spin Green’s func¬ 
tions D[l — l',t — t') and D^{1 — l',t — t') have been defined to describe properly 
spin propagations (SSl. 


3.2. Equation of motion method 


Since the spin operators obey the Pauli algebra, then our goal is to evaluate the 
charge-carrier and spin Green’s functions directly for the fermion and spin opera¬ 
tors in terms of the equation of motion method. In the framework of the equation 


of motion, the time-Fourier transform of the two-time Green’s function, 

G(u;) = ((A;At)), (40) 

satisfies the equation mg. 

w((yl; >!')>„ = {[A, A*]) + (([A, if]; At»„. (41) 

If the orthogonal operator L is defined as, [A,H] = — iL with ([L, = 0, 

the full Green’s function can be expressed as, 

G{u) = G(°)(cu) + 4 g(°)(u;)((L; L^))u.G{u;), (42) 


where ? = ([^djAl]), and the ME Green’s function, 

G(°)(u;) = 


uj-C 


(43) 


It has been shown 127 that if the self-energ y T ,(uj) is identified as the irreducible 
part of ((L;L'^))tj, the full Green’s function (42) can be evaluated as. 


G{u;) = 


with the self-energy. 


Ld — C — E(a;) ’ 


1 


S(a;) = -((L;Lt))-. 


(44) 


(45) 


In the framework of the diagrammatic technique, this self-energy E((u) corresponds 
to the contribution of irreducible diagrams. 


3.3. Mean-held approximation 


In the ME approximation, the t-J model (|34|) can be decoupled as |88| 
Hmfa 


Ht = 


Ht-\- - 2NZt(t)iXi + 2NZt'4>2X2, 

Xlt Y, - X2t' Y ^\+rahla - ^Y 


Ifja 


Ira 


la 


= X 


+2srsf+,i 


li) 


f<t'2Y.(StSi+f + SiS+ 




It 


(46a) 

(46b) 


(46c) 


18 










where the anisotropic parameter e = l+ 2 t</>i/Jeg, the charge-carrier’s particle-hole 
parameters (f>i = and (j )2 = {h\^hiJ^ro-)■, the spin correlation functions 

Xi = and X 2 = Z is the number of the nearest-neighbor or 

next nearest-neighbor sites on a square lattice, and N is the number of lattice 
sites. 

Within the framework of the equation of motion method, it is easy to find the 


MF charge-carrier Green’s function as 88 

g(0)(k,a;) = 


1 


a; - ?k 

where the MF charge-carrier excitation spectrum, 

?k = ^’fxiTk - ■^i'x27k - 


(47) 


(48) 


with 7 k = ( 1 /-Z^) = (cosA:a;-|-cosfcj;)/ 2 , andx^. = ( 1 /-Z^) = coskxcosky. 

However, in the doped regime without an AFLRO, i.e., {Sf) = 0, the MF spin 
Green’s func tions can be discussed in terms of the Kondo-Yamaji decoupling 
scheme [128 , which is a stage one-step further than the Tyablikov’s decoupling 
. After a straightforward calculation, the MF spin Green’s functions 


scheme 129 


have been o 


stained explicitly as 


, w = 


, w = 


2 wk 

B.{k) 


1 


1 


w — cjk a; -|- Wk 


1 


2a;z(k) \uj — a;z(k) 
with the MF spin excitation spectra. 




A? 
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1 


w -I- a;z(k) 
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(49a) 

(49b) 


Ml “ - aXiTk (e - 7k) A2 - ^^aexi - aex? 7 k 


(1 - e7k)] + Ai a ( xl7k “ ^^2 j 7k + ^ 


2Z 

^3 - 
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c^z^(k) = 


^ae(C'3 - X27k) 

eXl ^eAi - ^axi - axi7k^ (1 - 7k) + AM3(1 - 7k) 

+ AiA2[aeG3(7k - 1) + a(X27k - eC'3)(l - 7k)], 
and the functions, 

•Bk = Ai[2xi(e7k - 1) -t-Xi(7k - e)] - A 2 ( 2 xl 7 k - X 2 ), 

-Bz(k) = exiAi(7k - 1) - X2A2(7k - 1): 


(50a) 


(50b) 


(51a) 

(51b) 


where Ai = 2ZJeff, A 2 = 4 Z(^ 2 ^^ Ai = aCi -|- (1 — a)j{2Z), A 2 = aCf + (1 — 
a)/(4Z), A 3 = aC 2 + (1 ~ the spin correlation functions xf = {SfSf_^_~), 

xi = Cl = fi/z^) ,,>. Cf = (I/Z^) ^2 = 

(l/Z") Cs = (l/Z) Ef(Si7i+f>. “"d Cl = (l/Z) Ef(Sf+,Sf+f). 

In order to satisfy the sum rule of the correlation function (S'j^Sf) = 1/2 in the 
case without AFLRO, the i mpor tant decoupling parameter a has been introduced 
in the above calculation 88,128 , which can be regarded as the vertex correction. 
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3.4. Kinetic-energy driven superconductivity 


In the following discussio ns, we 1^-91 show that given the MF solution of the t-J 
model (34) in subsection 3.3 anciBy including the fluctuation around it due to the 


interaction between charge carriers and spins directly from the kinetic energy, we 
can obtain a formalism for the charge-carrier pairing which can be used to compute 
Tc on the first-principles as can be done for conventional superconductors. In 
particular, this formalism also can give a description of physical quantities which 
are consistent with the rather severe set by experiments. The interaction between 
charge carriers and spins directly from the kinetic energy in the t-J model (34) 
is quite strong, which can induce the SC-state in the particle-particle channel by 
the exc hange of spin excitations in the higher power of the doping concentration 
To see this point clearly, the self-consistent equations that are satisfied 


89-91 


by the full charge-carrier diagonal and off-diagonal Gre en’s funct ions are obtained 
in terms of the Eliashberg’s strong coupling formalism [130 131 


as. 


-h ff^°^(k,a;)[S^’"^(k,a;)c/(k,a;) - S^’"^(-k,-w)r'^(k, w)], (52a) 

rt(k,, 


= 


(-k, -w)[s[*"^(-k, -a;)r'^(-k, -w) -a;) 5 ((k,a;)], (52b) 


(h)/ 


where the self-energies S^*^^(k, w) in the particle-hole channel and Si^^^kjO;) in the 


(h). 


particle-particle channel are evaluated from the spin bubble as 89-91 


^P+P'-4-k « Z] 5r(p k, ipra + fa;n)n(p, p', ipra) 


_/\r2 “p+p'-r-K p 

P,P' 

I/eff(k,p,ZPm)5'(P + k,Z 


ipm + Wn 


(53a) 


iPm 


1 1 

S2^^(k,fa;n) = jVm) 


J^2 ^“P+P'+K^ 

p,p' ^Pm 

1 


N 


X] 4 X] ^eff(k, P, ipm)T\-T£> - k, 
P ipm 


'^Pm ) 1 


(53b) 


respectively, with the effective charge carrier interaction, 

I4ff(k,p,a;) = ^^A2_^p,+kn(p,p',a;), (54) 


and the spin bubble, 

n(p,p',ipm) = ^'^D^^Hp,ip'm)D^^Kp+Pjp'm + iPm)- 


(55) 


Wn 


This spin-excitation-mediated interaction (54) is a key to the SC transition in 
cuprate superconductors. Since both the pairing force and charge-carrier pair 

order parameter have been incorporated into the self-energy S 2 ^^(k,C(j), it is called 
the charge-carrier pair gap in the charge-carrier excitation spectrum, 


Ah(k,u;) = s('^)(k,a;), 


(56) 


which corresponds to the energy for breaking a charge-carrier pair and creating two 
charge-carrier quasiparticles. On the other hand, the self-energy S^*^^(k,a;) renor¬ 
malizes the MF charge-carrier spectrum, and therefore it describes the charge- 
carrier quasiparticle coherence. In particular, the self-energy ^^'^^(kja;) is an 
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even function of uj, while the self-energy w) is not. In this case, the self¬ 
energy can be broken up into its symmetric and antisymmetric parts as, 

S^^^(k, w) = S^g^(k, w) -|-a;E^Q^(k,a;), then both E^g^(k,a;) and S^Q^(k,a;) are an 
even function of ui. Moreover, the antisymmetric part w) of the self-energy 

S^^^(k, cj) is directly related to the charge-carrier quasiparticle coherent weight as, 

1 . _ ^(h). 


ZhF(k,a;) 


= l-ReSi“^(k,a;). 


(57) 


As a first step of discussions, we [^-91 study the kinetic-energy driven su¬ 
perconductivity, and therefore only focus on the low-energy behavior. In this 
case, the charge-carrier pair gap and quasiparticle coherent weight can be gen¬ 
erally discussed in the static limit, i.e., Ah(k) = S 2 ^^(k,t(j) |a;= 0 j and = 

I — ReS[Q^(k, w) |aj=o- As in conventional superconductor 

function ReS^^^ (k, cj) 1^=0 just renormalizes the chemical potential. Although 
ZhF(k) still is a function of momentum, the momentum dependence may be unim¬ 
portant in a qualitative discussion, and therefore the wave vector k in ZhF(k) can 
be chosen as. 


130 


131 


hF ' 

the retarded 


1 


=l-ReSS;;)(k,o; = 0) |k=[, 


, 0 ]) 


(58) 


133|. Moreover, this charge-carrier 


just as it has been done in the experiments 132 
quasiparticle coherent weight Zff reduces the charge-carrier (then electron) quasi¬ 
particle bandwidth, and suppresses the spec tral weight of the single-particle exci¬ 
tation spectrum, then the energy scale 121 of the quasiparticle band is controlled 
by the magnetic exchange coupling J. In particular, this charge-carrier quasi¬ 
particle coherence antagonizes superconductivity, and then is depressed to low 
temperatures 89 -91 . On the other hand, the s-wave component of the charge- 


carrier pair gap is suppressed heavily by Zff, and then the SC-state is dominated 
by the d-wave component. 


Ah(k) = Ah7k 


(d) 


(59) 


with 7 ^'^^ = (coskx—cosky) / 2 , which is consistent with the experimental fact 
that the charge-carrier pair state in cuprate superconductors has a dominant d 
wave symmetry over a wide range of the doping concentration, around the optimal 
doping. With the above static limit approximation, the full charge-carrier diagonal 
and off-diagonal Green’s functions in Eq. (52) are obtained explicitly as, 
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\a; - Ehk 

rt(k,u;) = 
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'^hk 
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2Ehk “ -®hk + -£^hk 


(60a) 

(60b) 


where Ehk = y I Ahz(k) P is the charge-carrier quasiparticle energy spec¬ 
trum, ^k = ■^hF'^k is the renormalized charge-carrier excitation spectrum, and 
Ahz(k) = ZhFAh(k) is the renormalized charge-carrier pair gap, while the charge- 
carrier quasiparticle coherence factors. 
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= F 1 1 - 


gk 

-Eihk 
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with the constraint = 1 for any wave vector k (normalization). In spite 

of the pairing mechanism driven by the kinetic energy by the exchange of spin ex¬ 
citations, the results in Eqs. (60) and (61) are the standard BCS expressions for a 
d-wave charge-carrier pair state. However, as a natural consequence of the charge- 
spin recombination, this charge-carrier pair state also leads to form the electron 
pairing state, and then the obtained d-wave BCS-like formalism for the electron 
pairing 89-91,119,134 indicates clearly the Bogoliubov quasiparticle nature of 
the SC quasiparticle peak. As in conventional superconductors, the Bogoliubov 
quasiparticle in cuprate superconductors does not carry definite charge, and is a 
coherent combination of the particle and its absence, then the SC coherence of 
low-energy excitations and the related quantities can be discussed on the first- 
principles basis much as can be done for conventional superconductors. Moreover, 
the AFSRO correlation has been incorporated into the SC-state through the spin’s 
order parameters entering into the charge-carrier self-energies (53) in the particle- 
particle and particle-hole channels, therefore there is a coexistence of the SC-state 
and AFSRO correlation, and then AFSRO fluctuation persists into superconduc¬ 
tivity. 


3.5. Self-consistent equations 
With the help of the full charge-carrier Green’s functions in Eq. 


Green’s function in Eq. (49a), the self-energies S^*^^(k, w) and ^^"''(k, w) are eval 
uated explicitly as (89 - 9 
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respectively, with u = 1,2, Bp+pU{^^ p'^p+p')^ ^/^pp' = 

<^p+p' ~ {—^Y^p'-, and the functions, 

-^iVpp'k ~ ?T-F(£'hp+k){l + ^T'BCa.'p'+p) + raB[(~l) ^^p']} 

+ nB(u;p/+p)reB[(-l)'''^+p'], (63a) 

■^ihip'k = [l-raF(-Ehp+k)]{l + nB(wp/+p)+nB[(-l)''++p']} 

+ nB(wp/+p)nB[(-l)'''^+p'], (63b) 

where nB{u}) and nF(ca) are the boson and fermion distribution functions, re¬ 
spectively. In this case, the charge-carrier quasiparticle coherent weight Zpp and 
charge-carrier pair gap parameter Ah satisfy following two self-consistent equa- 
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respectively, where kA = [tt, 0]. These two equations (64a) and ( 64b|) must be 
solved simultaneously with following self-consistent equations [89 -92,134 -136], 


4>i = 

4>2 = 

<5 = 

Xi = 

X2 = 

Cl = 
C2 = 

Cs = 

1 

2 


xl 


xl 


2N 


2N 


7k^hF f 1 

k ^ 

X] 7k-^hF (1 


^tanh[^/3Fhk] 

^hk ^ 

tanh[-/?Fhk] 


^hk 


k 

-V 

N 

k 

1 / ^k 

k 

vE- 


-Etanh[i/3i?hk] 

^hk ^ 


2^_^coth[i/3a,k|, 


k 

vEl'l'^cothli/Jcakl. 

1 

^2^7k7k 


N 


2 a,k”“'l 2 '®‘^l’ 


^E^”‘*'lsiS“kl, 


N ^ 2u;k 

k 


^ Z]^k;^^^coth[J/3a;^(k)], 


iV 


‘ 2 w 2 (k) 


./ ^.(k) 


^ IZ7k77^^coth[7/3a;^(k)], 


N 


‘ 2 w 2 (k) 


rtz — _ 

- N 




Cl = — 
3 iV 


-V 

V 


^ 2 a;,(k) 


(65a) 

(65b) 
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(65g) 

(65h) 

(65i) 

(65j) 

(65k) 

(651) 
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then all the order parameters, the decoupling parameter a, and the chemical po¬ 
tential /i are determined b y the self- consistent calculation without using any ad¬ 
justable parameters (^-92 ,134- 136[ . 

3.6. Doping dependence of charge-carrier pair gap and coupling strength 


Th e above e quations (64) and (65) have been calculated self-consistently , 90 


92,134-136 , and the result 135 of the charge-carrier pair gap parameter Ah as 


(64a) 

(64b) 
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Figure 2: The charge-carrier pair gap parameter as a function of doping with T = 0.002 J 
for t/J = 2.5 and t'/t = 0.3. [From Ref. 


135 


a function of doping for parameters t/J = 2.5 and t'/t = 0.3 with temperature 
T = 0.002J is shown in Fig. where the charge-carrier pair gap parameter 
Ah takes a domelike shape with the underdoped and overdoped regimes on each 
side of the optimal doping (^optimal ~ 0.15, where Ah reaches its maximum. The 
Andreev reflection experiments measure directly the binding energy of the charge- 
carrier pair while the Raman scattering and INS are the reverse experiments 
of the Andreev reflection in which they excite a charge-carrier pair out of the 
condensate energy [7, 16,29 . All the experimental data from the Andreev reflec¬ 
tion, Raman scattering and INS (ZllEI 2^[^ indicate that the charge-carrier pair 
gap parameter follows Tc as a function of doping. The domelike shape of the 
doping dependence of Ah in Fig. [^obtained within the framework of the kinetic- 
energy driven SC mechanism is qualitatively consistent with these experimental 

mi an 


28,29|. In particular, since the charge-carrier pair order is established 


data 

through an emergence of the charge-carrier quasiparticle, the charge-carrier pair 
state (then the SC-state) is controlled by both the charge-carrier pair gap Ah(k) 
and charge-carrier quasiparti cle c ohere nce ^ hF; which is reflected directly from 
the self-consistent equations (64a) and (64b). Moreover, this charge-carrier pair 
gap parameter Ah is strongly temperature dependent. To show this poin t clea rly, 
the charge-carrier pair gap parameter Ah as a function of temperature 136 for 
t/J = 2.5 and t'/t = 0.3 at doping 6 = 0.09 is shown in Fig. For co mpari- 
son, the corresponding experimental result of the pair gap parameter [137| for the 
underdoped Bi 2 Sr 2 Ca 2 Cu 30 io +5 is also shown in Fig. (inset). This calculated 
result in Fig. indicates that the charge-carrier pair gap parameter Ah follows 
qualitatively a BCS-type temperature dependence, i.e., it decreases with increasing 
temperatures, and eventuall y van ishes at Tc, which is also in qualitative agreement 
with the experimental data [137 of cuprate superconductors. 


In the framework of the kinetic-energy driven SC mechanism, spin excitations 
are directly coupled to charge-carrier pairs, and then the strength of this coupling 
with each charge-carrier pair is measured by the charge-carrier pair gap parameter 
2Ah. If the strength of the coupling exceeds the pairing energy 2Ah, the charge- 
carrier pairs will immediately be broken up. Since both the coupling strength Ves 
and charge-carrier pair order parameter Ah have been incorporated into the charge- 
carrier pair gap parameter Ah, the strength I4ff of the charge-carrier attractive 
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T/Tc 


Figure 3: The charge-carrier pair gap parameter as a function of temperature at ^ = 0.09 
for t/J = 2.5 and t'/t = 0.3. Inset: the corresponding experimental data of 
the pair gap parameter for the underdoped Bi2Sr2Ca2Cu30io+5 taken from 
[From Ref. 


Ref. 137 


136 


interaction mediated b y spin ex citations can be therefore obtained in terms of the 
ratio of Ah and A 


h as 


92 135 


VeS = 


( 66 ) 


where the charge-carrier pair order parameter A h can be evaluat ed explic itly from 
the charge-carrier off-diagonal Green’s function (60b) as 89-92,134-136 , 






(67) 


In Fig. we 92,135 show (a) the coupling strength Vgs and (b) charge-carrier 


pair order parameter Ah as a function of doping for t/J = 2.5 and t'/t = 0.3 with 
T = 0.002J. For co mparison, the correspond ing e xperimental results of the cou¬ 
pling strength _[138] and pair order parameter [139| for Bi2Sr2Ca2Cu308+5 are also 
shown in Fig. and Fig. El (inset), respectively, where the pair order parameter 
is related to the SC peak ratio (SPR), which is defined as the SC peak intensity di¬ 
vided by the overall spectral weight in the antinodal point, and then the pair orde r 
parameter is obtained indirectly by the experimental measurements of SPR [139 


The result in Fig. & shows that the coupling strength Ves smoothly decreases 
upon the increase of doping from a strong-coupling case in the underdoped regime 
to a weak-coupling side in the overdoped regime, and therefore is qu alitatively 
consistent with the experimental result of cuprate superconductors 12,138 140 . 
However, the charge-carrier pair order parameter A^ increases with Increasing 
doping in the lower doped regime, and reaches a maximum around the critical 
doping (5criticai ~ 0.195, then decreases with increasing doping in the higher doped 
regime. In comparison with the corresponding result of the charge-carrier pair gap 
parameter Ah in Fig. it is therefore found that the special doping dependence 
of the coupling strength I4fr in Fig- Eh induces an important shift from the critical 
doping ^critical ~ 0.195 for the maximal Ah to the optimal doping ^optimal ~ 0.15 
for the maximal Ah. 


25 




































5 



6 


Figure 4: (a) The coupling strength and (b) the charge-carrier pair order parameter as 
a function of doping with T = 0.002J for t/J = 2.5 and t'/t — 0.3. In¬ 
set in (a): the corresponding experimental data of the coupling strength for 

and inset in (b): the corresponding 


Bi2Sr2Ca2Cu308+5 taken from Ref. 138 


experimental data of the pair order parameter for Bi 2 Sr 2 Ca 2 Cu 30 g _|_5 taken 
[From Ref. 


from Ref. 139 


135 


3.7. Doping dependence of Tc 


Tc on the other hand can be obtained self-consistently from the self-c onsistent 
equations (64) and (65) by the condition Ah = 0, and the result |l35| of Tc as 
a function of doping for tjJ = 2.5, t'/t = 0.3, and J = lOOOK is shown in Fig. 
^ in comparison with the corresponding experimental results (inset). Since 
the charge-carrier pair state (then the SC-state) of cuprate superconductors is 
controlled by both the charge-carrier pair gap and charge-carrier quasiparticle 
coherence, which leads to that the maximal Tc occurs around the optimal doping, 
and then decreases in both the underdoped and the overdoped regime s, in good 
agreement with the experimental results of cuprate superconductors [l3 141| 142 


This calculated result also shows that the charge-carrier pairs (5^ are formed 


in momentum space, and then as in conventional superconductors 126 |130 131 


the pairing and onset of the phase coherence take place simultaneously at Tc- In 
particular, Tc that is set by the pair gap and quasi particle co herence has been 
observed experimentally in cuprate superconductors 132 133 , which may be a 


common feature for all superconductors. This follows from a fact that in spite 
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Figure 5: Tc as a function of doping for t/J = 2.5, t'/t = 0.3, and J = lOOOK. Inset: the 
corresponding experimental result of cuprate superconductors taken from 
[From Ref. 


135 


13 


in conventional s upercond uctors is also 
[130 , 131 . Finally, it 


of the electron-phonon SC mechanism, 
determined by the pair gap and quasiparticle coherence 
should be emphasized that except for the quasiparticle coherence, the superfluid 
density ps in cuprate superconductors is closely related to the charge-carrier pair 
order parameter A^, which is why ps increases with increasing doping in the lower 
doped regime, and reaches a maximum arou nd the critica l doping ^critical ~ 0 -195, 
then decreases in the higher doped regime 135[ 143[ 144 . This is also why [139 
the charge-carrier pair order parameter Ah can be obtained indirectly in terms of 
the experimental measurement of the superfluid density ps- We shall come back to 
give an interpretation about this issue in the discussions of the doping dependence 
of the superfluid density (see subsection 4.5). 

In the fermion-spin theory (23), the electron is decoupled as the charge carrier 
and spin according to its charge and spin degrees of freedom. However, we work 
in the case t 3> J, so that charge carrier motions are much faster than spins. 
Although the charge carriers repel each other because of the Coulomb interaction, 
at low energies there is an effective attraction resulting directly from the interaction 
between the charge carriers and spins in the kinetic energy of the t-J model ( |34[ ) 
by the exchange of spin excitations. This follows from a fact that a doped Mott 
insulator is formed by mobile charge carriers detaching themselves from spins 
that form the spin configuration. Such a mobile charge carrier is a magnetic 
dressing, and therefore its motion is strongly dependent on the spin configuration. 
When this charge carrier moves, the spin configuration rearrangements due to the 
presence of this charge carrier itself is left in its wake. This attracts a second 
charge carrier, leading to a net attraction between charge carriers. This kinetic- 
energy driven SC mechanism works because the spin dynamics is slow relative 
to charge carriers, which also is a natural consequence of the fact that the spin 

. In comparison with 
the spin excitation in 
e phonon in conventional 


excitations are very local in cuprate superconductors [24 
the conventional electron-phonon SC mechanism 34, 
cuprate superconductors plays a similar role to that of 
superconductors. 

On the other hand, the essential physics of a domelike shape of the doping 
dependence of Tc in the framework of the kinetic-energy driven SC mechanism can 
be attributed to a competition between the kinetic energy and magnetic energy in 
cuprate superconductors. The parent compounds of the cuprate superconductors 
are Mott insulators When charge carriers are doped into a Mott insulator. 
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there is a gain in the kinetic energy per charge carrier proportional to t due to 
hopping, however, at the same time, the magnetic energy decreases, costing an 
energy of approximately J per site [^. This leads to that the spin excitation 
spectral intensity decreases with increasing doping. However, a decrease of the 
spin excitation spectral intensity with increasing doping also leads to a decrease of 
the coupling strength Ves with increasing doping as shown in Fig. To see this 
competition clearly, (solid line), 5/ (5 m (dotted line), and 2Ah (dashed 

line) as a function of doping with T = 0.002J for tlJ = 2.5 and t'/t = 0.3 is 
shown in Fig. where = 14fr|5Ri0.045 is the value of Heff at the starting point 
of the SC dome, while (5max ~ 0.27 is the doping concentration at the end point 
of the SC dome. In the underdoped regime, the coupling strength Ves is very 
strong to bind the most charge carriers into charge-carrier pairs, and therefore 
the number of charge-carrier pairs increases with increasing doping, which leads 
to that the charge-carrier pair gap parameter Ah and Tc increase with increasing 
doping. However, in the overdoped regime, 14ff is relatively weak. In this case, 
not all charge carriers can be bound to form charge-carrier pairs by this weakly 
attractive interaction, and therefore the number of charge-carrier pairs decreases 
with increasing doping, which leads to that the charge-carrier pair gap parameter 
Ah and Tc decrease with increasing doping. In other words, in analogy to the 
electron-phonon SC mechanism for conventional superconductors, the reduction of 
Tc of cuprate superconductors in the overdoped side is driven by a reduction in the 
coupling strength I4fr of the pairing interaction. In particular, the optimal doping 
is a balanced point, where the number of charge-carrier pairs and coupling strength 
Ves are optimally matched. This is why the maximal Ah and Tc occur around the 
optimal doping, and then decrease in both the underdoped and overdoped regimes. 
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Figure 6: The coupling strength (solid line), the doping concentration (dotted line), and 
the charge-carrier pair gap parameter (2Ah) (dashed line) as a function of 
doping with T — 0.002J for tjJ = 2.5 and t'jt = 0.3. [From Ref. 
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3.8. Summary 

We have reviewed briefly the theory of the kinetic-energy driven superconductiv¬ 
ity. The kinetic-energy driven SC mechanism is purely electronic without phonons, 
where cuprate superconductors involve charge-carrier pairs bound together by the 
exchange of spin excitations. In particular, this kinetic-energy driven SC-state is 
conventional BCS-like with the d-wave symmetry, and then the obtained d-wave 
BCS type formalism for the charge-carrier pairing can be used to compute Tc and 
the related SC coherence of the low-energy excitations in cuprate superconductors 
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on the first-principles basis much as can be done for conventional superconductors, 
although the pairing mechanism is driven by the kinetic energy by the exchange of 
spin excitations in the higher powers of the doping concentration, and other exotic 
magnetic scattering [7,24 is beyond the d-wave BCS type formalism. Further¬ 
more, this kinetic-energy driven SC-state is controlled by both the charge-carrier 
pair gap and quasiparticle coherence, which leads to that the maximal Tc occurs 
around the optimal doping, and then decreases in both underdoped and overdoped 
regimes. The theory of kinetic-energy driven superconductivity also indicates that 
the strong electron correlation favors superconductivity, since the main ingredient 
is identihed into a charge-carrier pairing mechanism not from the external degree 
of freedom such as the phonon but rather solely from the internal spin degree of 
freedom of the electron. 


4. Electromagnetic response 


Superconductivity is characterized by exactly zero electrical resistance and ex¬ 
pulsion of magnetic fields occurring in superconductors when cooled below Tc. 
The later phenomenon is so-called Meissner effect |34], i.e., a superconductor is 


placed in an external magnetic field B smaller than tne upper critical held Be, the 
magnetic held B penetrates only to a penetration depth A (few hundred nm for 
cuprate superconductors at zero temperature) and is excluded from the main body 
of the system. This magnetic held penetration depth is a fundamental parameter 
of superconductors, and provides a rather direct measurement of the superhuid 
density ps {ps = 25 34 , which is proportional to the squared amplitude of 

the macroscopic wave function. In particular, the variation of the magnetic held 
penetration depth (then the superhuid density) as a function of doping and tem¬ 
perature gives the information about the nature of quasiparticle excitations and 
their dynamics |25| . Moreover, the magnetic held penetration depth can be also 
used as a probe or the pairing symmetry, since it can distinguish between a fully 
gapped and a nodal quasiparticle excitation spectrum 25,3^. The former results 


in the thermally activated (exponential) temperature dependence of the magnetic 
held penetration depth, whereas the latter one implies a power law behavior. This 
is why the hrst evidence of the d-wave pairing state in cuprate superconductors 
was obtained from the earlier experimental measurement of the magnetic held 


penetration depth 42 


Experimentally, by virtue of systematic studies using the muon-spin-rotation 
measurement technique, the essential feature of the evolution of the magnetic 
held penetration depth and superhuid density in cuprate superconductors with 
doping and temperature has been establish ed now for all the temperature T <Tc 
throughout the SC dome 42, 144 - 151|: (a) the magnetic held screening is 
found to be of exponential character |145 , T46] , in support of a local (London-type) 
nature of the electrodynamic response; (b) the magnetic held penetration depth is 
a linear temperature dependence at low temperatures except for the extremely low 
temperatures wher e a strong deviation from the linear characteristics (a nonlinear 
effect) appears [^,147,148 , which leads to that the superfluid density ps shows 
a crossover from the linear temperature dependence at low temperatures to a 
non linear one at the extremely low tem peratures as illustrated in Fig. 

150]; (c) the experimental measurements |144, 151 throughout the SC dome show 


that the superhuid density ps appears from the starting point of the SC dome, 
and then increases with increasing doping in the lower doped regime. However, 
this Ps reaches its highest value around the critical doping ^critical ~ 0.19, and 
then decreases at the higher doped regime, eventually disappearing together with 
superconductivity at the end of the SC dome (see Fig. |^. This is different from the 
case of the doping dependence of Tc, since the maximal Tc occurs around the the 
optimal doping (^optimal ~ 0.15, and then decreases in both the underdoped and the 
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8 12 
Temperature (K) 

Figure 7: (Color) The temperature dependence of the superfluid density measured for 
YBa2Cu306+5 at different doping concentrations. [From Ref. 


149 



Hole doping (5) 


Figure 8: The doping dependence of the low-temperature depolarization rate for un¬ 
derdoped to overdoped Tli_yPbySr 2 Cai_a;Ya;Cu 207 . The solid (dotted) line 
marks critical (optimal) doping. [From Ref. 


144 


overdoped regimes 
that the maximal 


13,141,142 


In particular, it has been shown (13 141,144 151| 
around the optimal doping (^optimal ~ 0.15 and the peak of 
Ps around the critical doping ^critical ~ 0.19 are a common feature of the hole- 
doped cuprate superconductors. Thus both Tc and ps variation as a function of 
doping provides important information crucial to understanding the details of the 
SC-state . Theoretically, the magnetic field penetration depth and the related 


superfluid oensity in cuprate superconductors have been discussed widely based 
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on a phenomenological d-wave BCS formalism [152 -156 . In the local limit, where 
the magnetic field penetration depth A is much larger than the coherence length 
C, i.e., A C) it ha-s been shown 153 that a simple d-wave pairing state gives the 
linear temperature dependence of the magnetic field penetration depth AA(T) = 
A(r) — A(0) oc r/Ao at low temperatures, where Aq is the zero-temperature value 
of the d-wave gap amplitude. However, the characteristic feature of the d-wave 
energy gap is the existence of the gap nodes, which can lead to the nonlinear 
effect of field on t he penetra tion depth (then superfluid density) at the extremely 
low temperatures [152 -156 . This follows from a fact that the nonlocal effect is 
closely related to the divergence of the coherence length ( at the gap nodes, as the 
coherence length varies in inverse proportion to the value of the energy gap, and 
then this nonlocal effect at the extremely low temperatures can lead to a nonlinear 
temperature dependence of the magnetic field penetration depth in the clean limit. 

In this section, we show how the kinetic-energy driven SC mechanism can 
be applie d to the discussions of the Meissner effect in cuprate superconductors 
[143 , |157 . For the discussions of the doping and temperature dependence of the 
electromagnetic response, the t-J model (34) can be extended by including the 


exponential Peierls factors as [135, 143,157 


H = 




l+f)/ 


li) 




It 

la If] 


( 68 ) 


where the exponential Peierls factors ac count for the coupling of the electron charge 
to an external magnetic field 158,159 in terms of the vector potential A(/). 


4.1. Linear response approach 

In a superconductor, an external magnetic field generally represents a large pertur¬ 
bation on a superconductor, then the induced field arising from the superconductor 
cancels this external magnetic field over most of the system. In this case, the net 
field acts only near the surface on a scale of the magnetic field penetration depth, 
and then it can be treated as a weak perturbation on the system as a whole . 
This is why the Meissner effect can be successfully studied within the linear re¬ 
sponse approach^! 


160 


161 


where the linear response current density Jn and the 


vector potential Ay are closely related by a kernel of the response function Kny as 


J^(q, uj) = -Y] K^y{ci, u;)Ay{ci, u), 


i/=i 


(69) 


where the Greek indices label the axes of the Cartesian coordinate system. It 
should be noted that the relation (69), which is local in the reciprocal space, in 
general implies a nonlocal response in the coordinate space. The kernel of the 
response function plays a crucial role for the discussion of the electromagnetic re¬ 


sponse, and can be separated into two parts as iL^jy(q, cj) = {q,io)+K^^ {q^,Lo), 


-(p)/ 


where and Kj^y^ are the corresponding diamagnetic and paramagnetic parts, 
respectively, and are related to the electron current density in the presence of 
the vector potential Ay. The evaluation of the diamagnetic contribution usually 
poses no difficulties since it is known almost immediately from the form of the 
diamagnetic current operator: it turns out to be diagonal and proportional to the 


Ap) 
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average kinetic term. However, the paramagnetic part can only be calculated ap¬ 
proximately since it involves evaluation of a retarded current-current correlation 
function (polarization bubble). 

The vector potential A (then the external magnetic held B = rot A) has been 
coupled to the electrons, which are now represented by Ci^ = hj^Sf and 6*4 = 
in the fermion-spin representation (23). However, in the CSS framework, the 

vector potential A is coupled to the electron charge in the t-J model ( 68 ), while 
the corresponding electron magnetic momentum can be coupled to the external 
magnetic held by including the Zeeman term in the t-J model ( | 68 [ ). However, for 
cuprate superconductors, t he upper critical m agnetic h eld is 50Tesla or greater 
around the optimal doping [162 163 , while we [143 157] mainly focus on the case 
where the applied external magnetic held (B < 10 mT) is much less than the 
upper critical magnetic held, and then the Zeeman term in the t-J model ( 68 ) has 
been dropped. We shall come back to give some comments about this issue in 
subsection |4.7| In this case, the electron current operator is obtained in terms of 


the electron polarization operator, which is a summation over all the particles and 
their positions, and can be expressed explicitly in the fermion-spin representation 
as [T4^[T^ , 


P = -e ^ KiClCi^ = e ^ Rih]hi, 

la I 


(70) 


then the electron current operator [131] is obtained by evaluating the time-derivative 
of the polarization operator as j = dP/dt = {i/h)[H, P]. In particular, in the lin¬ 
ear response approach, this electron current operator is reduced as j = 
with the corresponding diamagnetic (d) an d parama gnetic (p) components of the 
electron current operator that are given by 143, 157j, 


j(d) = 


eH 


^ f,A{i). f,{hi^hY^s^sr_^- + hi^hi-^sfsY) 


eH' 


If] 


— fA(0 . + hi^hl.^S^SY), (71a) 


?(p) - 


iet 


It 




l+fj’' 


Ifj 


iet' 

h 


Y+t)^ 


(71b) 


It 


respe ctively. Since the diamagnetic component of the electron current operator 
(71a) is proportional to the vector potential, the diamagnetic part of the response 
kernel is obtained directly as. 


- - 2x2</>2t')V = ('^2) 

where ~ /h^ is the doping and temperature dependence 

of the London penetration depth. 

However, the paramagnetic part of the response kernel is more complicated 
to calculate, since it involves evaluation of the following electron current-current 
correlation function. 


= -(rUiW(q,T)i<I')(-q,0)}), (73) 
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then the paramagnetic part of the response kernel (q, w) is obtained as (q, ui) 


cj). In the fermion-spin approach (23), the paramagnetic component of the 


electron current operator (|71b) can be decoupled as 


J"' = 


iex2t' 


h 

iecpit 


lf}(T 


It(j 


ie(j)2t' 


E + srsE) + ^ E ^(Esr+f+sr-six).(74) 


ifj 


It 


where the third and fourth terms in the right-hand side refer to the contribution 
from the electron spin, and can be expressed explicitly as. 


E »[(Esr+a +-sr-sE) - (Esr-s + srstj] 

ld=x,y 

= E “[(Esr+a + srsE)-EE'Sr+-sr+aE)i = t). psa) 

ld=x,y 

^ Ei(^+^(Esr+i+i+EE+x+i) - (i+!/)(Esr-i-i+-srEx-j) 

i 

+ a- !i)(s+s,-+ s,-s+ -(i- s)(s,+sr-i+«+srEi+pi 
= ^ Ek^*+ s)(EE+i+s+EE+i+s) - (i+s)(E+i+iE+E+*+sE) 

i 

+ (i- sXEE+i-i + EE+i-i) -a- !/)(E+i-«E + E+i-sE)l = ousb) 


which shows that the majority contribution for the paramagnetic component of the 
electron current operator comes from the electron charge, however the strong inter¬ 
play between charge carriers and spins has been considered through the spin’s order 
parameters entering in the charge-carrier part of the contribution to the current- 
current correlation. In particular, if we want to keep the theory gauge invariant, 
it is crucial to approxi mate the c orrelation function in a way of maintaining local 
charge conservation 34,159 -161 . Since the calculatio ns in this subsection will be 
worked on with a fixed gauge of the vector potential 143,157 , we postpone the 


detailed discussions of the gauge invariant problem until subsection |4.6[ Starting 


with the paramagnetic current operator (74), we can obtain its Fourier trans¬ 
form in the Nambu representation in terms of the charge-carrier Nambu operators 
'^ic ~ (^kf '^'k+q = (hk+qt; fel k-qj, )'^' purpose of addressing 

the gauge invariance problem in subsection |4.6| it is convenient to find the charge- 
carrier Green’s functions and electron density in the Nambu notation as well. From 
Eq. (60), the charge-carrier BCS-type Green’s function with the d-wave symmetry 
can be expressed in the Nambu representation as. 


g(k, io) = ZhF 


wro -b ^krs - Ahz(k)ri 


- ^^k 


(76) 


where tq is the unit matrix, ri and ra are the Pauli matrices. Since the density 
operator is summed over the position of all particles, its Fourier transform can be 

obtained as p{q) = (e/2) X)ka ^L^k+qa = (e/2) X)k ^k'^3^k+q, then the param¬ 
agnetic density-current operator can be represented in the Nambu representation 
as. 


J?Hq) = ^l7/.(k, k + q)4'k+q. (77) 

kcr 
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with the bare current vertex. 


7^(k + q,k) = < 


^ ~T e2*'^''{sin(A:^ + \q^,)[xit - 2X2^' E cos{\qy) cos{ku + \qu)] 

-i{2x2t') cos{k^ + 29m) Z) sin q^y sm{ku + \qv)]TQ for /r 7^ 0, 


i/7^M 


ers 


for /i = 0. 


It is necessary to be aware that we [143|1^ are calculating the polarization bubble 
with the paramagnetic current operator (|77[>, i.e., the bare current vertex (78), but 


the full charge-carrier Green’s function ( |76p . Consequently, as in th is scenario we 
do not take into account longitudinal excitations properly 34,159 , the obtained 
results are valid only in the gauge, where the vector potential is purely transverse, 
e.g. in the Coulomb gauge. In this case, the correlation function (73) can be 
obtained in the Nambu representation as, 

vE 7/,(k + q,k)7*(k + q, Tr[5(k-hq,'ia;„-hir'm)ff(k,ir'm)]. (79) 


Substituting the charge-carrier Green’s function (76) into Eq. ( |79p , the paramag¬ 
netic part of the response kernel in the static limit (w ~ 0) is obtained as. 


+ + +^2(k,q)] = i7jp)(q,0)<5/,^, (80) 


N 


where the functions Ei(k, q) and E2(k, q) are given by, 

72 / 1 I 6c+qCk + Ahz(k -I- q)Ahz(k)\ nF(Ehk) — ?^F(Ehk+q) 


-^i(k, q) — ^hF 1 1 + 
Z2(k, q) = ZIt? ( 1 - 


EhkEhk+q / Ehk — Ehk+q 

6c+q'?k + Ahz(k -I- q)A hz(k) ^ 2’'(Ehk) + ?^F(Ehk+q) “ 1 


EhkEhk+q 


Ehk + Ehk+q 


respectively. Now the kernel of the response function is obtained from Eqs. (72) 
and (80) as, 




i+JClP>(q. 0 ) 


lAI 




(82) 


(81a) 

(81b) 


4.2. Doping dependence of Meissner effect in long wavelength limit 

In the long wavelength limit, i.e., |q| —)• 0, the function L 2 (k, q —)• 0) vanishes, 
an d the n the paramagnetic part of the response kernel can be obtained explicitly 
as |143| , 

-f^yy^(q^ 0,0) = 2ZgF^^ ^ sin^ ky[xit - 2 x 2 t' cos k^f lim 


At zero temperature T = 0, it is found that Kyy^(q —)• 0, 0)|7’=o = 0, and then the 
long wavelength electromagnetic response is determined by the diamagnetic part 
of the kernel only. On the other hand, at T = Tc, the charge-carrier gap parameter 
Ah jT='Z} . = 0, and then the paramagnetic part of the response kernel is evaluated 
as |143| , 

-f^yy^(q ^ 0, 0)|T=rc = 2^hF^^ hbcit - 2 x 2 t' cos k^f lim — ^^(^k+q) 

N ^ q^O - 4k+q 


1 

'4’ 


(84) 


(78) 
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which exactly cancels the diamagnetic part of the response kernel (72), and then 


the Meissner effect in cuprate superconductors is obtained for all the temperatures 
T < Tc- To show this point clearly, the effective superfluid density ns{T) at 
temperature T is defined in terms of the paramagnetic part of the response kernel 


as, 


K(P)(q^0,0) = -^ 


Air 


1 - 


ns(r) 

ns( 0 ) 




and then the kernel of the response function ( |82| ) can be rewritten as, 

1 n (T) 


(85) 


( 86 ) 


where the ratio UsiT)/ns{0) of the effective superfluid densities at temperature T 
and zero temperature is given by. 


ns(r) 

ns(0) 


= 1 - 2A 


2 72 1 


^sin^ ky[xit 
k 


2 x 2 t' cos /Cj 


(e/3®hk -|- 1)2 


(87) 


In Fig. we show the effective superfluid density ns(r)/ns(0) as a function 
of temperature for t/J = 2.5, t'/t = 0.3, and J = lOOOK at 6 = 0.09 (solid 
line), 6 = 0.12 (dashed line), and 6 = 0.15 (dash-dotted line), where the effective 
superfluid density diminishes with increasing temperatures, and disappears at Tc, 
then all the charge carriers are in the normal fluid for the temperatures T >Tc. 



T(K) 


Figure 9: The effective superfluid density as a function of temperature at d = 0.09 (solid 
line), S = 0.12 (dashed line), and <5 = 0.15 (dash-dotted line) for t/J = 2.5, 
t'/t = 0.3, and J = lOOOK. [From Ref. 


143 


The main results obtained within the kinetic-energy driven SC mechanism 


are summarized as [l43l : (a) the doping dependence of the Meissner effect in 
cuprate superconductors is obtained for all the temperatures T < Tc throughout 
the SC dome; (b) the electromagnetic response kernel goes to the London form in 
the long wavelength limit [see, e.g., Eq. (86)]; (c) although the electromagnetic 
response kernel is not manifestly gauge invariant within the bare current vertex 
(|78|), however, we can keep the gauge invariance with in the dressed current vertex 
157 , which will be proven clearly in subsection 4.6 
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4.3. Quantitative characteristics 


The way the system reacts to an external electromagnetic stimulus is entirely 
described by the linear response kernel (69), which is calculated within the kinetic- 
energy driven SC mechanism. Once the response kernel is known, the effect 
of an electromagnetic field can be quantitatively characterized by experimentally 
measurable quantities such as the local magnetic field profile and the magnetic 
field penetration depth. Technically, it is needed to combine one of the Maxwell 
equations with the relation (69) describing the response of the system and solve 
them together for the vector potential. This is the step in which a particular 
gauge of the vector potential — usually implied by the geometry of the system — 
is set. However, the result we have obtained the response kernel (82) can not be 
used for a direct comparison with the corresponding experimental data of cuprate 
superconductors because the kernel function derived within the linear response 
theory describes the response of an infinite system. In order to take into account 
the confined geometry of cuprate superconductors, it is necessary to introduce 
a surface being the boundary between the environment a nd the sa mple. This 
can be done within the standard specular reflection model 164,165 with a two- 
dimensional geometry of the SC plane, in the configuration with external magnetic 
field perpendicular to the ab plane, as illustrated in Fig. [T^ In this subsection we 
study magnetic field penetration effects within the ab plane only, so our goal is to 
find and discuss the in-plane magnetic field penetration depth. 



X - penetration depth 


Figure 10: The geometry of the specular reflection model. The current Jext simulates 
external magnetic field at the edge of the sample (a; = 0), whereas the induced 
supercurrent Jint is the (linear) reaction of the system. [From Ref. 


157 


In order to simulate an external magnetic field at the surface of a two-dimensional 
sample, we introduce an external current sheet Jy,ext(a^) = —2B(5(x)//ro at the edge 
X = 0, where /tq is the magnetic permeability and B is the amplitude of the exter¬ 
nal magnetic field B at the surface (x = 0). From the Maxwell equation for the 
curl of the local magnetic field roth = /tq (Jint+Jext) = /^oJint + [0, —2B(5(x), 0] and 
the fact, that the induced supercurrent Jint flows along the y axis, we can state 
that the local magnetic field is of the form h(r) = [0, 0, hfix)]. In order to discuss 
the magnetic field penetration effect, the spatial dependence of the local magnetic 
field has to be found. Let us begin with the identity rot rot A = grad div A — A 
and choose the vector potential as A(r) = [0, Ay(x), 0] setting the Coulomb gauge. 
In this case, q^Ay{(fi) = /io[Jy,int(q) + Jy,ext(q)]j because the vector potential has 
only non-zero y component. Finally, including the form of the external current, 
the linear relation (69) between the induced supercurrent and the vector potential 
Jy,int(q) = ~Ayy(q)Ay(q), and solving for the vector potential we obtain. 


^y(q) = -2B 


^oAyy(q) ql' 


( 88 ) 
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Since the vector potential has only the y component, the only non-zero component 
of the local magnetic field h = rot A is that along the z axis and /iz(q) = f<?xAy(q). 
Substituting the derived form of the vector potential (88), and taking the inverse 
Fourier transform, the local magnetic held prohle can oe obtained explicitly as, 


K{x) = — 

TT 


dqx 


Qx sm{qxx) 


f^O^yyiqxj 0 , 0 ) + q^ 


(89) 


This local magnetic held pro hles can be measured experimentally, e.g. using the 
/iSR measurement technique 145,146 , providing an important tool to investigate 


the details of magnetic held screening inside the sample. For convenience in the 
following discussions at this section, we introduce a characteristic length scale 
oo = . Using the lattice parameter a « 0.383nm for YBa 2 Cu 307 _ 2 y, 

this characteristic length is obtain as oq ~ 97.8nm. In this case, the local magnetic 
held prohle (89) as a func tion of th e distance from the surface has been studied 
at different doping levels [143 , 157 , and the theoretical results perfectly follow 
an exponential held decay as expected for the local electrodynamic response. In 
particular, this exponential character of the local magnetic held prohle has been 


observed experimentally on different families of cuprate superconductors 145 146 
in support of a local (London-type) nature of the electrodynamics. 


4.4. Doping and temperature dependence of magnetic held penetration depth 


The local magnetic held prohle hz{x) in the Meissner state obtained in Eq. (89) 
allows us to determine the magnetic held penetration depth A(r) in a straightfor- 
ward way. According to the dehnition A(T) = hz{x) dx, the magnetic held 

penetration depth can be evaluated as. 


MT) = 1 


hz{x)dx Tj; J ^^Kyy{qx,0,0) + q^' 

0 


(90) 


At zero temperature, the calculated magnetic held penetration depths 143 157 
are A(0) ss 239.17nm, A(0) ps 234.76nm, and A(0) « 224.44nm at <5 = 0.14, 
5 = 0.15, and <5 = 0.18, respectively, which are qualitatively consistent with the 
values of the magnetic held penetration depth A ~ 156nm ~ 400nm observed 


for different fami lies of cuprate superconductors at different doping levels [l44 


145,149, 166, 167 . On the other hand, at T = Tc, the kernel of the response 


function Km)(q ^ 0,0)|r=rc = 0, and then the magnetic held penetration depth 


from Eq. (90) can be found as A(rc) = oo, which rehects that in the normal- 
state, the external magnetic held can penetrate through the main body of the 
syst em, therefore there is no the Meissner effect in the normal-state. In Fig. [H 
we 111 show the magnetic held penetration depth AA(T) = A(T) — A(0) as a 
function of temperature for t/J = 2.5, t'/t = 0.3, and J = lOOOK at <5 = 0.14 
(solid line), d = 0.15 (dashed line), and 5 = 0. 18 (d ash-dotted line) in comparison 


with the corresponding experimental results 

(H 3 of YBa 2 Cu 307 _y (inset). In 

low temperatures, the magnetic held penetration depth AX{T) exhibits a linear 
temperature dependence, however, it crosses over to a nonlinear behavior in the 
extremely low temperatures, in good agreement with expe rimental observation 
in nominally clean crystals of cuprate superconductors (4^ p5Hl47| . However, 


it should be emphasized that the result in Fig. 11 for cuprate superconductors 


is different from that in conventional superconductors, where the characteristic 
feature is the existence of the isotropic SC gap Ag, and then AA(T) exhibits an 
exponential behavior as AX{T) oc exp(—Ag/T). 
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Figure 11: The temperature dependence of the magnetic field penetration depth AA(T) 
at 5 = 0.14 (solid line), ^ = 0.15 (dashed line), and 6 = 0.18 (dash-dotted 
line) for t/J = 2.5, t'/t = 0.3, and J = lOOOK. Inset: the corresponding 

p^ . [From Ref. [I43 


experimental data for YBa 2 Cu 307 _y taken from Ref. 


4.5. Doping and temperature dependence of superfluid density 

The superfluid density Ps(T) can be obtained directly from the magnetic field 
penetration depth A(T) as, 


PsiT) = 


A2(r)' 


(91) 


The zero temperature superfluid density PsjO) as a function of doping for t/ J = 2.5, 
t'/t = 0.3, and J = lO OOK is shown in Fig. 12 in comparison with the correspond¬ 
ing experimental data 1144 for Yo, 8 Cao, 2 Ba 2 pui_^Zn 2 ) 307_5 and Tli_j^Pbj^Sr 2 Cai_ 
(inset). This calculated result clearly shows that the superfluid density Ps{0) in¬ 
creases with increasing doping in the lower doped regime, and reaches a high¬ 
est value (a peak) around the critical doping 6 ~ 0.195, then decreases in the 
higher doped regime. In particular, this anticipated value of the critical doping 
^critical ~ 0.195 is Very close to the critical doping (^critical ^ 0-19 observed e xperi- 
mentally for different families of cuprate superconductors [144 , 149, 151[ 166 . The 
early experimental data observed from cuprate superconductors show that the su¬ 
perfluid density /9s(0) in the underdoped regime vanis hes m ore or less linearly with 
decrease of the charge-carrier doping concentration 5 (nS). Later, a clear deviation 
from this linear relation between the superfluid density /Os(0) an d charge- carrier 
doping concentration has been observed in th e und erdoped regime |144|148 . How¬ 
ever, the recent experimental measurement [149 on the cuprate superconductor 
YBa 2 Cu 307 _y indicate that the superfluid density /9s(0) is, in actual fact, lin¬ 
early proportional to the charge-carrier doping concentration in the doped range 
5 ~ 0.054 ~ 0.061. The calculated result in Fig. a in this doped range also is 


149 


well consistent with the experimental observation 

However, the critical doping ^critical ~ 0.195 for the highest ps is different 
from that in the doping dependence of Tc as shown in Fig. and the doping 


xYxCu207 
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Figure 12: The doping dependence of the superfluid density in T = 0 for t/J = 
2.5, t'/t = 0.3, and J = lOOOK. Inset: the corresponding ex¬ 
perimental results for Yo. 8 Cao. 2 Ba 2 (Cui_ 2 Zn 2 ) 307_5 (open circles) and 
Tli_yPbySr 2 Cai_a;Ya;Cu 207 (solid triangles) taken from Ref. [g. [From 


Ref. 
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dependence of the charge-carrier pair gap parameter as shown in Fig. 
where the maximal and Ah appear around the optimal doping ^optimal ~ 0.15. 
This difference is a long-standing puzzle in cuprate superconductors, however, its 
interpretati on is quit e simple within the framework of the kinetic-energy driven SC 
mechanism [135 143] . This follows from a fact that the superfluid density Ps from 
Eqs. (91) and (|90pls closely related to the kernel of the response function (82), and 
therefore the charge-carrier pair order parameter Ah, the coupling strength 14 ff, 
and all the other order parameters are relevant. In other words, the variation of 
the superfluid density ps with doping and temperature is coupled to the d oping and 
temperature dependence of Ah, I4g, and all the other order parameters 135,143 . 
In particular, the doping-derivative of ps at the critical doping (^critical ~ 0.195 
is obtained as (d/9s/d(5)|5=5^^j,;^^; = 0. Since p^ = A-^, (dps/d<5)|5=5^^;,;^^; = 0 is 
equivalent to (dA/d(5)|^=^ , = 0. In this case, (dA/d(5)|5=5^^j,.j^^j = 0 can be 

expressed in terms of Eq. (|90| as. 


'dA' 


/ 

dqx 

1 dKyy(qx,0,0)' 

.d5. 

^—^critical 

^ J 

[/ioKyy(qx,0,0)-kq|]2 dJ 


6=Sc: 


= 0 , 


and then it is straightforward to obtain from Eq. (82) that when (dpg/d(5) = 

0, (dAh/d(5)|5=5^^j,.j^^j = 0, which shows that the^ping effects from the coupling 
strength I4g and all the other order parameters upon p^ are almost canceled each 
other, and then the behavior of the doping dependence of the superfluid density 
Ps is mainly dominated by the doping dependence of the charge-carrier pair order 
parameter Ag. However, the charge-carrier pair order parameter Ag measures the 
strength of the binding of two charge carriers into a charge-carrier pair, and it has 
been shown in subsection |3.6| that Ag has a domelike shape of the doping depen¬ 
dence with the maximal value appearing around the critical doping ^critical ~ 0.195 
(see Eig. ^). In particular, the charge-carrier pair order parameter Ag and the 
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charge-carrier pair macroscopic wave functions in cuprate superconductors are the 
same as we have mentioned in section i.e., the charge-carrier pair order param¬ 
eter is a magnified version of the charge-carrier pair macroscopic wave functions. 
On the other hand, the superfluid density ps is a measurement of the phase stiff¬ 
ness [^, and is proportional to the squared amplitude of the charge-carrier pair 
macroscopic wave functions. Both ps and Ah thus describe the different aspects of 
the same charge-carrier pair macroscopic wave functions. In this case, the domelike 
shape of the doping dependence of ps with the highest value appearing around the 
critical doping is a natural consequence of the domelike shape of the doping depen¬ 
dence of Ah with the maximal value appearing around the same critical doping. 
In comparison with the results obtained in subsection 3.7 it is therefore shown 


that except for the quasiparticle coherence, ps is determined by the charge-carrier 
pair order parameter Ah, while is set by the charge-carrier pair gap parameter 
Ah, this is why there is a difference between the optimal doping for the maximal 
Tc and the critical doping for the highest ps in cuprate superconductors. 



T(K) 


Figure 13: The temperature dependence of the superfluid density at 6 = 0.06 (solid line), 
S = 0.09 (dashed line), S = 0.12 (dash-dotted line), and 6 = 0.15 (dotted 
line) for t/J = 2.5, t'/t = 0.3, and J = lOOOK. Inset: the corresponding 
experimental result for YBa 2 Cu 307 _y taken from Ref. [149| .[From Ref. 


143 


The doping dependence of the superfluid density shown in Fig. 12 is also 


strongly temperature dependent. In particular, when the temperature T = Tc, 
the kernel of the response function itr^i/(q —)• 0,0)|r=rc = 0, and then the mag¬ 
netic field penetration depth A(rc) = oo as mentioned above, which leads to the 
superfluid density Ps{Tc) = 0, which is consistent with t he re sult of the effective 
superfluid density obtained from Eq. ( [ST] ). In Fig. 13, we [143] show the superfluid 
density pfiT) as a function of temperature for tj J = 2.5, t' jt = 0.3, and J = lOOOK 
at 5 = 0.06 (solid line), 5 = 0.09 (dashed line), <5 = 0.12 (dash-dotted line), and 
5 = 0.15 (dotted line) in comparison with the corresponding experimental re¬ 
sult H of YBa 2 Cu 307 _y (inset). pfiT) decreases with increasing temperature, 
and vanishes at Tc. In particular, the most striking feature of the results is the 
wide range of linear temperature dependence at low temperatures, extending from 
close to Tc to down to the temperatures T ~ 4K~ 8K for different doping concen- 
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trations. However, in correspondence with the nonlinear temperature dependence 
of the magnetic field penetration depth at the extremely low temperatures shown 
in Fig. [m the superfluid density Ps{T) crosses over to a nonlinear temperature 
behavior at the extremely low temperatures (below T ~ 4K~ 8K for different dop¬ 
ing concentrations). These calculated re sults are also qualitatively consistent with 
the corresponding experimental results [144,149[ 166| of cuprate superconductors. 


The explanation [143 157 for the nonlinearity in the temperature dependence 
of the penetration depth (then the superfluid density) at the extremely low tem¬ 
peratures i s the sam e as the case based on the phenomenological d-wave BCS 
formalism [152 - 156] , and can be found from the nonlocal effects induced by the 
gap nodes m a pure d-wave pairing state. An external magnetic field acts on 
the SC-state of cuprate superconductors as a perturbation. Within the linear 
response theory, one can find that the nonlocal relation between the supercur¬ 
rent and the vector potential (69) in the coordinate space holds due to the finite 
size of charge-carrier pairs. In particular, in the kinetic-energy driven d-wave SC 
mechanism (M -91 , the size of charge-carrier pairs in the clean limit is of the or¬ 
der of the coherence length ^(k) = fiup/vrAh (k) , where up = is 

the charge-carrier velocity, which shows that the size of charge-carrier pairs is 
momentum dependent. In general, although the external magnetic field decays 
exponentially on the scale of the magnetic field penetration length A(T), any non¬ 
local contributions to measurable quantities are of the order of where the 
Ginzburg-Landau parameter k is the ratio of the magnetic field penetration depth 
A and the coherence length (. However, for cuprate superconductors, because the 
pairing is d-wave as shown in Eq. (59), the charge-carrier pair gap vanishes on 
the gap nodes, so that the quasipartiHe excitations are gapless and therefore af¬ 
fect particularly the physical properties at the extremely low temperatures. This 
gapless quasiparticle excitation leads to a divergence of the coherence length ("(k) 
around the gap nodes, and then the behavior of the temperature dependence of 
the magnetic field penetration depth (then the superfluid density) depends sen¬ 
sitively on the quasiparticle scattering. At the extremely low temperatures, the 
quasiparticles selectively locate around the gap nodal region, and then the ma¬ 
jor contribution to measurable quantities comes from these quasiparticles. In this 
case, the Ginzburg-Landau ratio ^(k) around the gap nodal region is no longer 
large enough for the system to belong to the cl ass o f type-H superconductors, and 
the condition of the local limit is not satisfied [T^ , which leads to the system in 
the extreme nonlocal limit, and therefore the nonlinear behavior in the tempera¬ 
ture dependence of the magnet ic field p enetration depth (then superfluid density) 
is observed experimentally 25,145,146[. On the other hand, with increasing tem¬ 


peratures, the quasiparticles around the gap nodal region become excited out of 
the condensate, and then the nonlocal effect fades away, where the momentum 
dependent coherence length C (k) can be replaced approximately with the isotropic 
one Co = ^up/vrAh. In this case, the calculated Ginzburg-Landau parameters are 
kq « A(0)/Co ~ 166.29, kq « 175.55, and kq « 156.14 for the doping concentra¬ 
tions (5 = 0.14, (5 = 0.15, and 6 = 0.18, respectively, and then the condition for 
the local limit is satisfied. In particular, these calculated values of the Ginzburg- 
Landau parameter at different doping concentrations are very close to the range 
Ko ~ 150 ~ 400 estimated experiment ally for different famil ies of cuprate super¬ 
conductors at di&rent doping levels [144, 145,149 


166 167 . As a consequence, 


the study 143,157 based on the kinetic-energy driven SC mechanism shows that 


cuprate superconductors at moderately low temperatures turn out to be type-H 
superconductors, where nonlocal effects can be neglected, and then the electrody¬ 
namics is purely local and the magnetic field decays exponentially over a length 
of the order of a few hundreds nm. 
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4.6. Gauge-invariant electromagnetic response 


Although the electromagnetic response kernel is not manifestly gauge invariant 
within the bare current vertex (78), however, we can keep the theory gauge in- 

[T^ . 


variance within the dressed current vertex 


It is 


invariance is a direct consequence of local charge conservation 34,160 


well known that gauge 
which is 


mathematically expressed by the charge density-current contin uity equatio n or its 
Green function analogue called the generalized Ward identity 


34 


159-161 


2 

+ Q,k) = T3g~^{k) - g-^{k + q)T^, 

fi=0 


(93) 


where the charge e and reduced Planck constant h have been set to the unity, P^ is 
a dressed version of the density-current vertex function, and for convenience in the 
discussions at this subsection, the three-vector notation q = (q, qo = iujn) along 
with the metric (1,1, —1) has been introduced. 

Since the local charge conservation requirement is quite universal and funda¬ 
mental, it should be inherent to any theory of the Meissner effect which is expected 
to be gauge invariant. The purpose of this subsection is to propose-within the for¬ 
malism of kinetic-energy driven superconductivity-a method to dress the current 
vertex in a way, which does not violate the generalized Ward identity. Once such a 
method is found, the bare polarization bubble (73) can be replaced by its dressed 
version presented in Fig. 14, and the resulting kernel of the response function will 


provide correct results for any gauge of the vector potential. 


Y 



9 


r 


Figure 14: Dressed polarization bubble (Nambu notation). Here both the charge-carrier 
Green function and the current vertex are dressed with the pairing interaction 
due to the spin bubble. [From Ref. [157 .] 

In the first step we will note that [157 
2 

'^q^ll^^{k + q,k) = T^g^^^-^{k) - g^^'i-^ik + q)T2., (94) 

/r=0 


i.e. the generalized Ward identity for the bare current vertex is satisfied with 
the MF charge-carrier Green function (47) g^^\k) = [{iujnf' — {ioJnTo + 
Substituting this MF charge-carrier Green function, the right-hand side of Eq. (94) 
turns into T 3 g^^^~^{k) — g^^^~^{k + q)T^ = (Ck+q ~ Ck)'^o ~ Moreover, in the 
long wavelength limih after including the explicit form of the MF charge-carrier 
dispersion relation (48), it further simplifies to T 3 g^^^~^{k) — g^^^~^{k + q)T^ ~ 
[-qxsm.kx{ 2 txi - Z¥x2 cosky) - Qy sin/c^(2tyi - Zt'x2COskx)]To - q^T^. Now, 
recalling the form of the bare vertex (78), it is easy to find that in the long 
wavelength limit the scalar product on the left-hand side of Eq. (94) is equal 
to the above obtained one of the right-hand side of Eq. (94), which proves the 
equality (94). 
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Figure 15: Ladder-type approximation for the dressed vertex. [From Ref. 


157 .] 


It is well known that in order to obtain a dressed vertex function, which 
does not violate the generalized Ward identity, a ladder-type approximation can 


be adapted 34,159 ,1601. The nature of the pairing mechanism 89,90 , which 


form 157 


originates from the spin bubble (55), suggests a ladder-like approximation of the 


1 1 


r^{k + q,k) = -f,^{k + q,k) + —-J2'^39{k+p + q)Tf,{k+p + q,k + p)g{k + p)T3 


^ ^P+p'+kn(p, p', ipm), 


(95) 


which is graphically presented in Fig. 15 


In order to prove that the approximation (95) for the dressed vertex in fact 


implies a gauge invariant description of the Meissner effect, it is necessary and 
sufficient to check whether it does not violate the generalized Ward identity ( |93[ ). 
In order to prove it, we insert the dressed vertex function (95) into the left-hand 
side of Eq. ( [j^ and use the identity Yll^=o^^^^^^(k + q, k) = T‘ig~^{k)—g~^{k + q)T‘i 
to obtain 


1 1 




p=0 


p=0 

1 


Y ^p+p'+kn(p, P , iPm). 


(96) 


p' 


In the long wavelength limit we use the approximation ss Ap^p,^j^^q for 

the second term of the right-hand side in Eq. (96). Then we can simplify Eq. 


(96) in terms of the self-energy (53) into + 9’^) 


_ _ ^ _ Ep=o9/^7p(fe + 

q, k) + + q)T 3 — T^'E^^^k). Using the fact that the bare vertex satisfies the 

generalized Ward identity with the ME charge-carrier Green function, as st ated 
in Eq. ( |94[ ), and arranging the terms with respect to the Pauli matrices, we [157 
have 


YQ^^^,ik + q,k) ^ Ts[g^^^-\k)-E(-^\k)] - [gio)-\k + q) - E^^^k + q)]rU97) 


^ 1=0 


Hence, identifying the terms in the square brackets as the full charge-carrier Green 
functions, we eventually obtain the generalized Ward identity (93), which proves 


that the ladder-type approximation (95) for the vertex function in the dressed 
polarization bubble in Fig. 14 is consistent with the generalized Ward identity. 


Gonsequently, the kernel of the linear response calculated with the dressed polar¬ 
ization bubble is gauge invariant. 
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4.7. Summary and discussions 


Within the framework of the kinetic-energy driven SC mechanism, the doping 
dependence of the electromagnetic response is discussed. In the linear response 
approach, the electromagnetic response consists of two parts, the diamagnetic cur¬ 
rent, which is the acceleration in the magnetic field, and the paramagnetic current, 
which is a perturbation response of the excited quasiparticle and exactly cancels 
out the diamagnetic term in the normal state, then the Meissner effect is obtained 
for all the temperatures T < Tc throughout the SC dome. By considering the 
two-dimensional geometry of cuprate superconductors within the specular reflec¬ 
tion model, the main features of the doping dependence of the local magnetic 
field profile, the magnetic field penetration depth, and the superfluid density are 
qualitatively reproduced. The local magnetic field profile follows an exponential 
law, while the magnetic field penetration depth shows a crossover from the linear 
temperature dependence at low temperatures to a nonlinear one at the extremely 
low temperatures. In particular, the domelike shape of the doping dependence of 
the superfluid density ps with the highest value appearing around the critical dop¬ 
ing (^critical ~ 0.195 is a natural consequence of the domelike shape of the doping 
dependence of the charge-carrier pair order parameter Ah with the maximal value 
appearing around the same critical doping. 



Figure 16: The superfluid density as a function of temperature at 5 = 0.09 with magnetic 
field B = 0 (solid line), i? = 0.5 T (dashed line), and B = 1.0 T (dash-dotted 
line) for tjJ = 2.5, t'jt — 0.3, and J = lOOOK. Insets: the corresponding 
experimental results for YBa2Cu3 06.95 taken from Ref. |170| . [From Ref. 

T^.] 


Finally, it should be emphasized again that in the above discussions, the only 
coupling of the electron charge to an external magnetic field is considered in the 
t-J model (68) in terms of the vector potential A, while the coupling of the elec¬ 
tron magnetic momentum with the external magnetic field in terms of the Zee- 
man mechanism has been dropped. However, the depairing due to the Pauli spin 
polarization is very important in the presence of a moderate or strong external 
magnetic field, since cuprate superconductors are doped Mott insulators with the 
strong AFSRO correlation dominating the entire SC phase [7,24 . In particular, 
within the framework of the kinetic-energy driven SC mechanism, a moderate or 
strong external magnetic field aligns the spins of the unpaired electrons, then the 
d-wave charge-carrier pairs in cuprate superconductors can not take advantage of 
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the lower energy offered by a spin-polarized state [168 . In this case, the magnetic 


field de pende nce of the superfluid density in cuprate superconductors has been 
studied [169 by considering both couplings of the electron charge and electron 
magnetic momentum with a weak magnetic field, and the calculated result of the 
superfluid density as a function of temperature at <5 = 0.09 with the magnetic 
field i? = 0 (solid line), B = 0.5 T (dashed line), and B = 1.0 T (dash-dotted 
line) for t/J = 2.5, t'/t = 0.3, and J = lOOOK is shown in Fig. 


16 


m comparison 


with the corresponding experimental data El for YBa 2 Cu 3 06.95 (inset). Most 
importantly, the magnitude of ps at the extremely low temperatures decreases with 
increasing magnetic field, and then it turns to be independent on a weak magnetic 
field away from the extremely low temperatur es, in qua litative agreement with 
experimental data of cuprate superconductors 170- 173j. The calculated result 


also indicates that the nature of the quasiparticle excitations at the extremely low 
temperatures is strongly influenced by a weak magnetic field. This weak magnetic 
field induced reduction of the superfluid density of cuprate superconductors at 
the extremely lo w tem peratures contrasts with that observed from conventional 
superconductors 174 , where the curves of the temperature dependent superfluid 


density for differently weak magnetic fields were found to collapse onto a single 
curve since conventional superconductors are fully gaped. 


5. Dynamical spin response 

As illustrated in the schematic phase diagram in Fig. cuprate superconductors 
exist on a continuum with a family of the Mott insulators. The undoped and 
extremely low-doped cuprates have an AFLRO at low temperatures. The spin 
excitation with AFLRO is called as magnon. However, in the doped regime of the 
SC dome, the AFLRO breaks down, but the AFSRO correlation remains. This 
AFSRO can still support spin waves, but the spin excitations with AFSRO are 
damped. The damped spin excitation is known as paramagnon [125 . 


1.8 
1.6 

1.4 

1.2 

UJ 

0.8 

0.6 

0.4 

0.2 
0.0 

Figure 17: (Color) The magnetic dispersion relation in different families of cuprate su¬ 
perconductors. [From Ref. [^.] 
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Figure 18: (Color) The resonance energy as a function of doping referenced to the opti¬ 
mal doping, (5opt, corresponding to The red full curve shows the doping 

dependence of times 5.3. [From Ref. 


185 


The early INS measurements [7 23,175 - 183| on cuprate superconductors have 
demonstrated that the doped charge carriers cause substantial changes to the low- 
energy spin excitation spectrum, and a consistent pattern has been identihed as 
the hour-glass-shaped dispersion as illustrated in Fig. This hour-gla ss-sh aped 

and 


182 


dispersion was hrst obse rved in the spin excitations of YBa 2 Cu 306.6 
Lai.875Bao.i25Cu04 (ill. where two IC components of the low-energy spin excita¬ 
tion spectrum are separated by a commensurate resonance energy at the waist of 
the hour glass. In the upward component, above the commensurate resonance en¬ 
ergy (Ur, the spin excitation spectrum is similar to what one would expect from AF 
spin fluctuations with a finite gap, and is relevant to the results for different fam¬ 
ilies of cuprate superconductors that appear to scale with the magnetic exchange 
coupling constant J for the undoped parent compounds of cuprate superconduc¬ 
tors. In particular, for a given excitation energy, the magnetic scattering peaks 
lie on a circle of radius of with the incommensurability parameter that 
is dehned as a deviation of the peak position from the AF wave vector [1/2,1/2] 
(for convenience, in this section we use the units of [27r, 27r]) in the Brillouin zone 
(BZ), and then the distribution of the spectral weight of IC magnetic scattering 
peaks is rather isotropic. On the other hand, in the downward component, below 
ujr, the distributio n of the s pectral weight of IC magnetic scattering peaks is quite 
anisotropic 7, 23 175- 183|. In particular, it is remarkable [7 


23 


175 

the 


189 that 


in analogy to the domelike shape of the doping dependence m Tc, tne commen¬ 
surate resonance energy uj^ increases with increasing doping in the underdoped 
regime, and reaches a maximum around the optimal doping, then decreases in 
the overdoped regime as illustrated in Fig. 18, reflecting a intrinsical relationship 
between Wr and Tc. Although the IC magnetic scattering has been also observed 
in the normal-state, the c ommensu rate resonance is a new feature that appears 
in the SC-state only [7,23,175-189 . Later, this hour-glass-sh aped disp ersion was 
found in several different families of cuprate superconductors 


190 


193] . However, 

because of technical limitations, only the low-energy (T ~ 10 — 80 meV) spin 
excitations in a small range of mom entum sp ace around the AF wave vector are 
detected by INS measurements [23, 175- 193]. In recent years, instrumentation 


for RIXS with both soft and hard X-rays has improved dramatically, allowing this 
technique to directly measure the high-energy (T ~ 80 — 500 meV) spin excitations 
of cuprate superconductors in the w ide e nergy-momentum window that cannot be 
detected by INS measurements [^, 194 . In this case, as a compensation for the 
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miss of a significant part of the spectral weight of spin excitations in INS stud' 


les 


[^,175-193 , the RIXS technique has been used to measure the high-energy 


spin excitations of cu prate sup erconductors in the whole doping range, and the 
experimental data 24,195 -198| indicates that the key feature of the high-energy 


spin excitations even in the overdoped regime is strikingly similar to that of the 

- 201| (see Fig. 


undoped parent compounds 199 


19). 


x=0.00 x=0.11 x=0.16 x=0.26 x=0.40 



Figure 19: (Color) The dispersion of the magnetic excitations in La 2 -a;Sra;Cu 04 as func¬ 


tion of Q and x. [From Ref. 196 


Within the framework of the kinetic-energy driven SC mechanism, the dynam¬ 
ical spin respons e of cupra te superconductors from low-energy to high-energy has 
been studied [88,202-204 , where one of the main results is that both the damped 


but well-defined dispersive low-energy and high-energy spin excitations exist across 


the whole doping phase diagram. In the SC-state 202j , the low-energy spin excita¬ 
tions are strongly renormalized due to the interaction between charge carriers and 
spins to form an hour-glass-shaped dispersion. In particular, the commensurate 
resonance is closely related to the process of the creation of charge-carrier pairs, 
and appears in the SC-state only, while the low-energy IC magnetic scattering is 
mainly associated with mobile charge-carrier quasiparticles, and therefore persists 
into the normal-state. On the other hand, the ch arge- carrier doping has a more 
modest effect on the high-energy spin excitations 202] , and the high-energy spin 
fluctuations bear a s triking resemblance to those found in the undoped parent 
compounds |199 -201 . In this section, we summarize a f ew calculate d results for 


the dynamical spin response of cuprate superconductors 88,202-204 


47 



































































5.1. Dynamical spin structure factor in superconducting-state 


In the framework of the CSS fermion-spin theory (23), the scattering of spins 
due to the charge-carrier fluctuation is responsible to the spin dynamics. For the 
discussions of the dynamical spin response in cuprate superconductor s, it is needed 


to calculate the full spin Green’s function, which can be expressed as 88,202-204 

1 


D(k,a;) = 


D(o)-i(k,a;) - S(«)(k,a;)' 


(98) 


In the SC-state, the spin fluctuation occurs in the charge-carrier quasiparticle 
background, and then the spin self-energy can be obtained within the framework 
of the equation of motion method in terms of the colle ctive charge-carrier modes 
in the particle-hole and particle-particle channels |202| as. 


+ ^p-l-q-l-k) ^ E + iqm + ipm) 


pq 


'tQm 


X [nS(P> iQm) - n{f]^(p, q, iqm)], 


where the charge-carrier bubble Ilgg (p, q, iqm) in the particl e-hole channel is ob¬ 
tained from the full charge-carrier diagonal Green’s function (60a) as. 


(99) 


n^g(P,qC9m) = ■^'^9{p,i<^n)g{p + <l,iuin + 'iqm), 


( 100 ) 


and is closely related to mobile charge-carrier quasiparticles, while the charge- 
carrier bubble npp(p, q, ig^) in the particle-part icle channel is obtained from the 
full charge-carrier off-diagonal Green’s function (60b) as, 


npr(p,q, igm) = ^ y~]ri(p,iu;^)r(p q, iuin + iqm), 


( 101 ) 


and therefore is directly associated with the creation of charge-carrier pairs. Sub¬ 
stituting the full charge- carri er Gre en’s function (60) and the MF spin Green’s 
function (49a) int o Eq s. (100), (101), and (99), the spin self-energy can be evalu¬ 
ated explicitly as H, 


S(^)(k,.;) = -^^(-l)-+ia(^) 


^+(p,q)^ij(k,p,q) 


pq,y 


+ 


-^-(p,q)-p',;-(k,p,q 


kpq 

(s). 


- [Wq+k - (-l)‘'’''^(-E'hp+q - -Fhp)]^ 


~ [“^q-l-k ~ (~l)^'’'^(-£ihp+q + -®hp)]^ 


( 102 ) 


where p = 1,2, = ^hF(^Lp + ^p+q+k)-®q+k/(2wq+k), the charge-carrier 

coherence factors for the processes, 


-^-f(p,q) = 1 + 

-^-(p,q) = 1- 


Cp^p-fq - Ahz(p)Ahz(p + q) 

-FhpF'hp-fq 

gp^p+q - Ahz(p)Ahz(p + q) 

EhpEhp+q 


(103a) 

(103b) 
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and the functions 


-P’i?(k,p,q) 




[^q+k ( 1) (-^hp+q -^'hp)]{^B(^q+k) [^F(-^'hp) ^F(-®hp+q)] 

(-l)"+inF[(-l)"Shp]nF[(-l)"+'^hp+q]}, (104a) 

[l^q+k - (-l)‘^’'’^(-E'hp+q + -£^hp)]{?T'B(Wq+k)[l “ ?^F(-Ehp) “ ''^F(-E'hp+q)] 

(-l)"+inF[(-l)"+'Shp]nF[(-l)"+'i?hp+q]}. (104b) 


With the help of the full spin Green’s function ( |98[), t he dynamical spin structure 
factor of cuprate superconductors is obtained as l202|, 


S'(k, uj) 


—2[1 + nB(w)]ImZl(k, oj) 

2[l + nB(a;)]S2lmS(")(k,a;) 

“ [w2 -ujI- SkReSW (k, w)]2 + [5kImE(^) (k, t<;)]2 ’ 


(105) 


where ImS(®^(k, w) and Re S^’^^(k , oj) are the corresponding imaginary and real 
parts of the spin self-energy (102), respectively. 


5.2. Universal Low-energy spin excitation spectrum 



k 


Figure 20: The energy dependence of the position of magnetic scattering peaks at d = 
0.21 with T = 0.002J for t/J = 2.5 and t'/t = 0.3. [From Ref. ’ 


202 


First of all, we discuss unusual feature of the low-energy magnetic scattering. 
Of course, at half-filling, the undoped parent compounds of cuprate superconduc¬ 
tors are Mott insulators, and then AFLRO gives rise to a commensurate peak 
at [1/2,1/2]. However, the calculated dynamical spin structure factor spectrum 
shows that when AFLRO is suppressed with doping, two IC magnetic scatter¬ 
ing modes separated by a commensurate resonance energy Wr are developed [202 . 
Well above the magnetic resonance energy the IC magnetic scattering peaks 
lie uniformly on a circle of radius of and then the distribution of the spectral 
weight of IC magnetic scattering peaks is quite isotropic. However, the geom¬ 
etry of the magnetic scattering is energy dependent. In particular, below Wr, 
although some IC satellite peaks appear along the diagonal direction of BZ, the 
main weight of IC magnetic scattering peaks is in the parallel direction, which 
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leads to an rather anisotropic distribution of the spectral weight of IC magnetic 
scattering peaks below To show the energy dependence of the position of the 
low-energy magnet ic sca ttering peaks clearly, the evolution of magnetic scattering 
peaks with energy 202 at <5 = 0.21 with T = 0.002J for t/J = 2.5 and t'/t = 0.3 
is shown in Fig. 20, where the hour-glass-shaped dispersion of the low-energy 


magnetic scattering peaks observed fr om different families of cuprate supercon¬ 
ductors is qualitatively reproduced 


7 23 


175-193 . In particular, in contrast to 


the case at energies below Wr, the spin excitations at the energies above Wr disperse 
almost linearly with energy, which has been observed experimentally on cuprate 


superconductors [7,23, 175 


193 


5.3. Doping dependence of commensurate resonance 


we 


Th e com mensurate resonance energy Wr is strongly doping dependent. In Fig. 

202[ show the commensurate resonance energy uj^ as a function of doping with 


T = 0.002J for t/J = 2.5 and t'/t = 0.3, where in analogy to the domelike shape of 
the doping dependence of Tc, the maximal Wr occurs around the optimal doping, 
and then decreases in both the underdoped and the o verdoped regimes, also in 
qualitative agreement with the experimental results 7, |185 , [I^ . In particular, 
using a reasonably estimative value of J ~ 100 meV, the anticipated resonance 
energy Wr = 0.46J ~ 46 meV in the optimal doping Jopt = 0.15 is not too far from 
th e resonance energ y uJj- ss 41 meV observed in the optimally doped YBa2Cu306+5 


aiB 180||185[|186| . This commensurate resonan ce ha s been also discussed based 


on the slave-boson approach at zero temperature [205 . 
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Figure 21: The magnetic resonance energy uj^ as a function of doping with T = 0.002J 
for tjJ = 2.5 and t'jt = 0.3. [From Ref. 


202 


5.4. Evolution of high-energy spin excitations with doping 

For a comparison of the high-energy spin excitati ons at d ifferent doping levels 
just as it has been done in RIXS experiments (24[ 195f 198] , the dynamical spin 
structure factor 5(k, w) as a function of energy along the k = [0, 0] to k = [0.5, 0] 
direction of BZ with T = 0.002J for t/ J = 2.5 and t'/t = 0.3 at 5 = 0.04 (non-SC 
regime), 6 = 0.09, 6 = 0.15, <5 = 0.21, and 6 = 0.25 is shown in Fig. [221 These 


calculated results [202] capture the qualitative feature of the high- energy sp in 


ex¬ 


citations observed experimentally on cuprate superconductors [^,195-198 . The 


high-energy spin excitations persist across the whole doping phase diagram with 
comparable spectral weight and similar energies, i.e., in contrast to the dramatic 
change of the low-energy spin excitations with doping, the high-energy spin ex¬ 
citations retain roughly constant energy as a function of doping, and the shapes 


50 
























































5=0.04 

5=0.09 

5=0.15 

5=0.21 

5=0.25 

[0.20, 0] 

[0.20, 0] 

[0.20, 0] 

[0.20, 0] 

[0.20, 0] 

■ A 


'....TV. 



[0.30, 0] 


[0.30, 0] 

[0.30, 0] 

[0.30, 0] 

[0.30, 0] 

[0.35, 0] 


[0.35, 0] 

[0.35, 0] 

[0.35, 0] 

[0.35, 0] 

[0.40, 0] 


[0.40, 0] 

[0.40, 0] 

[0.40, 0] 

[0.40, 0] 

[0.50, 0] 


[0.50, 0] 

[0.50, 0] 

[0.50, 0] 

[0.50,0] 


i 




- /\ 


■ I . 1. 1,1 > -J . I . . ... . n —, I ■ I . I . I—rn ^ —r i ■ i . I i —T i ■ i 

012345012345012345012345012345 


CO U 


Figure 22: The dynamical spin structure factor S'(k, w) as a function of energy along 
the k = [0, 0] to k = [0.5,0] direction of the Brillouin zone at <5 = 0.04, 
5 = 0.09, (5 = 0.15, <5 = 0^, and S = 0.25 with T = 0.002J for t/J = 2.5 and 
t' jt = 0.3. [From Ref. 


202 


of the high-energy magnetic scattering peaks in the heavily overdoped regime are 
very similar to those in the extremely low-doped and underdoped regimes, al¬ 
though the width of the high-energy spin excitations increases continuously with 
doping, consistent with the spin excitation being damped by the increasing doping. 
Furthermore, for example, the magnetic scattering peak on energy scale of 4.7J 
appears in the k = [0.2,0] point at b = 0.21, while the peak on the energy scale of 
5.0 J emerges in the k = [0.3,0] point, reflecting the dispersive nature of the high- 
energy spin excitations along the k = [0,0] to k = [0.5, 0] direction. In particular, 
this dispersion relation of the high-energy spin excitations in the overdoped regime 
resembles those in the extremely low-doped and underdoped regimes. 


5.5. Dispersion of spin excitations 


To determine the overall spin excitation spectrum Es(k) in Eq. (105), spin exci¬ 
tations at different momenta in the whole doping phase diagram have been deter- 
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Figure 23: (a) The dispersion of the spin excitations along the high symmetry directions 
of the Brillouin zone at d = 0.21 for t/J = 2.5 and t'/J = 0.3 with T = 
0.002J. (b) The experimental result of the dispersion of the spin excitations 
along the high symmetry directions of the Brillouin zone for YBa2Cu306+5 
taken from Ref. 


194 . [From Ref. 202 


mined by a self-consistent calculation [202 , 

a;2 = w^ + .BkReS(®)(k,w), 


(106) 


and the result shows that spin excitations are well defined at all momenta. In 


Fig. [2^, we 202 show £'s(k) as a function of momentum along the high symme¬ 


try directions of BZ at d = 0.21 with T = 0.002J for t/J = 2.5 and t'/t = 0.3. 
For comparison, the experimental result |194 of Es(k) along the high symmetry 
directions of BZ for the cuprate superconductor YBa2Cu306+5 is shown in Fig. 
2^ . It is shown that these theoretical calculations reproduce qual itatively the 
overall dispersion of spin excitations in cuprate superconductors [194 . In compar¬ 


ison with the spin excitation spectrum (the spin wave) of the parent compounds 


of cuprate superconductors 01101 199| - [20l] , the spin excitation spectrum in the 
doped regime has been renormalized due to the presence of the interaction be¬ 
tween charge carriers and spins directly from the kinetic energy of the t-J model 
(34). However, the charge-carrier doping does not uniformly renormalizes the dis¬ 


persion of spin excitations. In particular, the low-energy magnetic correlation is 
strongly reorganized, where two IC components of the spin excitation spectrum 
are separated by the commensurate resonance energy cjr, which therefore leads 
to an hour-glass-shaped dispersion of magnetic scattering peaks 
shown in Fig 


20 


Moreover, the spin excitations at energies well above w, 


193 


as 


disperse 


almost linearly with energy, which is similar to spin wave with a finite gap, reflect¬ 
ing a fact that the charge carrier doping strongly renormalizes the spin excitations 
at energies below to,-, but has a modest effect on the spin excitation dispersion 
at energies above Wr [^. However, in contrast to the case of the low-energy spin 
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excitations, the dispersion of the high-energy spin excitations in cuprate super¬ 
conductors in t he overdoped regime is strikingly similar to that of their parent 
compounds 194,199- 


201 


The physical interpretation |202 to the above results can be found from the 


special property of the spin self-energy w) in Eq. (102) obtained directly 

from the interaction between charge carriers and spins in the kinetic energy of the 
t-J model (34). Th is follows from a fact that the dynamical spin structure factor 
S'(k, oj) in Eq. (105) has a well-defined resonance character, where S'(k, oj) exhibits 
peaks when the incoming neutron energy u) is equal to the spin excitation energy 
Es(k), i.e., 


2 2 
a; 


- 5k^ReS(^)(kc,u;) = ~ 0, 


(107) 


for certain critical wave vectors kc, the magnetic scattering peaks appear, and then 
the weights of these peaks are dominated by the inverse of the imaginary part of 
the spin self-energy l/ImS(®)(kc,a;). In other words, the positions of magnetic 
scattering peaks are determined by both the spin excitation energy ii's(k) and the 
imaginary part of the spin self-energy ImS(®)(kc, w). At half-filling, the low-energy 
magnetic scattering peak locates at the AE wave vector [1/2,1/2], so the commen¬ 
surate AF peak appears there. However, away from half-hlling, the doped charge 
carriers disturb the AF background. In particular, within the framework of the 
kinetic-energy driven SC mechanism, as a result of the self-consistent interplay 
between charge carriers and spin s, the unusual magnetic correlation is de velop ed. 


5.1 


the spin self-energy S(®)(k,w) in Eq. (102) is 


As mentioned in subsection 
obtained in terms of the full charge-carrier diagonal Green’s function (60a) and 


off-diagonal Green’s function (60b), and renormalizes spin excitations. However, 


in the charge-carrier quasiparticle spectrum Ehk = 




Ahz(k) P in the full 

charg e-carrier diagonal Green’s function (60a) and off-diagonal_ Green’s function 
(60b), the maximal pk| appears around the nodal regio n, and ^k has an effective 
band width ITh ~ 2 J at the end of the SG dome [202 . In particular, this band 
width decreases with decreasing doping. However, the d-wave charge-carrier pair 
gap Ah(k) vanishes on the gap nodes, while the d-wave charge-carrier pair gap 
parameter Ah has a domelike shape of the doping dependence with the maximal 
Ah ~ 0.2J appearing around the optimal doping. These properties of ^k and 
Ah(k) lead to that the effective band width of the charge-carrier quasipa rticle 
spectrum Ehk is almost the same as ^kj and then the spin self-energy in Eq. (102) 
strongly renormalizes the spin excitation at the energies below Hh) but has a weak 
effect on the spin excitation at the energies above ITh- This is why the magnetic 
correlation at the energies below Wh is strongly reorganized, while the high-energy 
spin fluctuation bears a striking resemblance to those fonnd in the nndoped parent 
compounds. Furthermore, as seen from the spin self-energy (102), there are two 
parts of charge-carrier quasiparticles contribution to the spin self-energy reno rmal- 
ization. The contribntion from the first term of the right-hand side in Eq. (102) 


mainly comes from mobile cha rge-ca rrier quasiparticles, and the coherence factor 
for this process is given in Eq. (103a). This process mainly leads to the low-energy 


IC magnetic scattering, which can persist into the normal-state. However, th e ad- 
ditional contribution from the second term of the right-hand side in Eq. (102) 
originates from the creation of charge-c arrier pairs, and the coherence factor for 
this additional process is given in Eq. (103b). This additional process occurs in 


the SC -stat e only, and gives a dominant contribution to the commensurate resO' 
nance (175 


reflecting that the commensurate resonance is intimately related to 
superconductivity, and then appe ars in the SC-state only. In particular, it is easy 
to find from Eqs. (105) and (102) when the incoming neutron energy uj ~ 2Ah at 


the AF wave vector [1/2,1/2J, the commensurate resonance peak appears, which 
leads to that show the same domelike shape of the doping dependence as Tc. 
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This universal relationship oj^ ~ 2Ah between the resonance energy and charge- 
carrier pair gap par ameter in cuprate sup erconductors has been confirmed by the 
experimental data (7l|175[|180[|185||186[|206[ . The theory of the kinetic-energy driven 
SC mechanism thus naturally explains why the commensurate resonance is seen 
below Tc only, and why its doping and temperature dependences scale with the 
SC order parameter. 


5.6. Dynamical spin response in normal-state 


At the temperature T > Tc, the charge-carrier pair gap parameter = 0, and 
superconductivity disappears, then the system is in the norm al-state. In this 
the dynamical spin structure factor 5(k, cj) in Eq. (105) in the SC-state 


case, 


is reduced to that in the normal-state, where the spin self-energy is obtained in 
terms of the collective charge-carrier mode i n the par ticle-hole channel only, and 
the calculated results are summarized as (M, 202 - '204| : (a) the commensurate res¬ 


onance that originates from the creation of charge-carrier pairs and appears in 


the SC-state 180, 185,1^ is absent from the normal-state, while only 

the low-energy IC spin fluctuation in the SC-state persists into the normal-state. 
Moreover, the low-energy IC magnetic scattering peaks lie on a circle of radius Jic- 
Although some IC satellite peaks along the diagonal direction appear, the main 
weight of IC magnetic scattering peaks is in the parallel direction as in the case of 
the SC-state. Since the height of IC magnetic scattering peaks is determined by 
damping, the IC magnetic scattering peaks broaden and weaken in amplitude as 
the energy increases. In particular, the dynamical spin structure factor spectrum 
has b^en used to extract the doping dependence of the incommensurability param¬ 
eter djC) and the results show clearly that die increases progressively with doping 
at the lower doped regime, bu t saturat es at the higher doped regime, in qualitative 
agreement with experiments 179, 207] ; (b) Although the high-energy spin excita¬ 
tions retain roughly constant energy as a function of doping, the width of these 
high-energy spin excitations increases with increasing doping. In particular, the 
high-energy spin excitations, in their overall dispersion, their spectral weight, and 
the shapes of the magnetic scattering peaks, are strikingly similar to those in the 
corresponding SC-state, although the magnetic scattering peak in the normal-state 
is softenin g and bro adening, also in qualitative agreement with the experimental 
results 
particu’ 


R 


195 - 198) ; (c) the integrated dynamical spin susceptibility appears to be 
y universal, and is scaled approximately as oc arctan[aia;/T +a 3 {uj/T)^ 
where ai and 0 3 are constants, which is qualitatively consistent with the experi' 
ments |208 ,|209| 


5.7. Summary 

In this section, we have summarized some calculated results of the dynamical 
spin response of cuprate superconductors obtained based on the kinetic-energy 
driven SC mechanism. The spin self-energy in the SC-state is evaluated explicitly 
in terms of the collective charge-carrier modes in the particle-hole and particle- 
particle channels, and then employed to calculate the dynamical spin structure 
factor. The calculated results show the existence of damped but well-defined dis¬ 
persive spin excitations in the whole doping phase diagram. In particular, the 
low-energy spin excitations in the SC-state have an hour-glass-shaped dispersion, 
with commensurate resonance that originates from the process of the creation of 
charge-carrier pairs, and appears in the SC-state only, while the low-energy IC 
spin fluctuation is dominated by the process from mobile charge-carrier quasi¬ 
particles, and therefore can persist into the normal-state. The high-energy spin 
excitations in the SC-state on the other hand retain roughly constant energy as a 
function of doping, with spectral weights and dispersion relations comparable to 
those in the corresponding normal-state, although the magnetic scattering peak in 
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the normal-state is softening and broadening. The dynamical spin response probes 
the local magnetic fluctuation and is a very detailed and strin gent test of the mi¬ 
croscopic theory of superconductivity. The calculated results 202f 204 based on 
the kinetic-energy driven SC mechanism lead to the behaviors similar to that seen 
in the experiments. 


6. Theory of normal-state pseudogap state 


The discovery of superconductivity in cuprate superconductors has been under¬ 
lined by two salient phenomena, their high temperature superconductivity and the 
occurrence of the normal-state pseudogap [l4] - [22j . In particular, the normal-state 
pseudogap is particularly obvious in the undeidoped regime, where the charge- 
carrier concentration is too low for the optimal superconductivity. A number 
of experimental probes [14 -22 show that below a characteristic temperature T*, 


which can be well above Tc in the underdoped regime, the physical response of 
cuprate superconductors can be interpreted in terms of the formation of a normal- 
state pseudogap by which it means a suppression of the spectral weight of the 
low-energy excitation spectrum. Moreover, this normal-state pseudogap crossover 
temperature T* decreases with increasing doping in the underdoped regime and 
since Tc rises with doping, then T* seems to merge with Tc in the overdoped 
regime, eventually disappearing together with superconductivity at the end of the 
SC dome [T4] - [T6] , which leads to an anomalous normal-state pseudogap state at 
the underdoped regime and eventually a crossover to the normal-metal phase at 
the heavily overdoped regime. 

In the early days of Mott insulators, Mott [210 introduced the term pseudogap 
to indicate a minimum in the density of states at the Fermi energy, resulting from 
the Coulomb repulsion between electrons at the same site. During the last two 
decades, several scenarios were proposed to explain the formation of the normal- 
state pse udogap in cuprate superconductors. In particular, it has been argued 
[211 , 212| that the normal-state pseudogap originates from preformed pairs at T *, 
which would then condense (that is, become phase coherent) at Tc. On the other 
hand, it has been suggested that the normal-state pseudogap is distinct from the 
SC gap and related with a certain order which competes with superconductivity 
(2B- |215|. Furthermore, it has been proposed that the normal-state pseudogap 


is a combination of a quantum disordered d-wave superconductor and an entirely 


different form of competing order, originating from the particle-hole channel [216 
However, up to now, no general consensus for the normal-state pseudogap has 
been reached yet on its origin, its role in the onset of superconductivity itself, and 
not eve n on its evolution across the phase diagram of cuprate superconductors 
In particular, the key questions surrounding the normal-state pseudogap 


16 


phenomenon and its relevance to superconductivity have been raised |15[|16] : (a) 
What phase diagram is the correct phase diagram with respect to the normal- 
state pseudogap line? (b) Is the normal-state pseudogap the result of some one- 
particle band structure effect? (c) Is there a true order parameter defining the 
existence of a normal-state pseudogap phase? (d) Do the normal-state pseudogap 
and SC gap coexist? (e) Is the normal-state pseudogap a necessary ingredient for 
superconductivity? 

In section]^ we (89-91 have discussed the kinetic-energy driven SC mecha¬ 


nism, and shown that the interaction between charge carriers and spins directly 
from the kinetic energy by the exchange of spin excitations in higher powers of the 
doping concentration generates SC-state in the parti cle-p article channel. Based 
on this kinetic-energy driven SC mechanism, we (^,217 have developed a mi¬ 
croscopic theory of the normal-state pseudogap state, and shown that the same 
charge-carrier interaction arising through the exchange of spin excitations that gen¬ 
erates the SC-state in the particle-particle channel also induces the normal-state 
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pseudogap state in the particle-hole channel, indicating that the spin excitation 
plays a decisive role in formation of both the SC-state and normal-state pseudogap 
state. In this section, we review briefly this microscopic theory of the normal-state 
pseudogap state. 


6.1. Relationship between normal-state pseudogap and quasiparticle coherence 


Within the framework of the kinetic-energy driven SC mechanism, the normal- 

state pseudogap opens due to the strong electron correlation without symmetry 

fh') 

breaking. This follows from the fact that the charge-carrier self-energy ^k, cj) 


in the particle-hole channel in Eq. 
as 1^, 


(62a) also can be rewritten approximately 


sS'^)(k,a;) 


[2Apg(k)]^ 

U) “h -/Vfk 


(108) 


where is the energy spectrum of (k, uj). As in the case of the charge-carrier 

pair gap mentioned in subsection |3.4[ the interaction force and order parameter in 
the characterization of the normal-state pseudogap state have been incorporated 
into Apg(k), 


then in this case, Apg(k) is so-called the normal-state pseudogap. It 


should be emphasized that the equation (108) is an identity only in the case of 
w = 0, however, it is a proper approximation for the case of a; / 0. Since the 
normal-state pseudogap s[*^^(k,a;) originates from the charge-carrier self-energy 

S^^^(k, (u) in the particle-hole channel, it can be identified as being a region of the 
self-energy effect in the particle-hole channel in which the normal-state pseudogap 
suppresses the spectral weight of the low-energy excitation spectrum. In particular, 
this normal-state pseudogap is directly related to the charge-carrier quasiparticle 
coherent weight (58) as. 


ZhF 


= 1 _ RpV^^Vk ( ) = ni I =1-1- I 

i (.K) W Oj |k=[7r,0] T j^2 lk=[7r,0]! 


(109) 


which shows that the main effect of the normal-state pseudogap has been contained 
in the quasipa rticle coherent weight. As a byproduct, it is therefore established 
a relationship (109) between the normal-state pseudogap Apg and charge-carrier 
quasiparticle coherent weight Zhp. Since the SC-state in the kinetic-energy driven 
SC mechanism is controlled by both the SC gap and quasiparticle coherence as 
mentioned in subsection |3.6[ in this sense, the normal-state pseudogap is a nec¬ 
essary ingredient for superconductivity. However, there is a strong competition 


between the charge-carrier quasiparticle coherence and charge-carrier pairs 90 
which leads in cuprate superconductors to be reduced to lower temperatures, 
indicating that the normal-state pseudogap has a competitive role in engendering 
superconductivity. 


6.2. Interplay between superconductivity and normal-state pseudogap state 


Substituting the self-energy s[*^^(k, w) in Eq. (108) into Eq. (52), the full charge' 


carrier diagonal and off-diagonal Green’s functions can be obtained by considering 
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the interplay between the SC gap and normal-state pseudogap as, 
g(k,u}) = 


w - Ck - - A 2 (k)/[a; -a;)] 


(h). 


Tj'i 

^Ihk 


-h 


T/2 

•^Ibk 


^ — -E'lhk ^ + -Elhk 


-h 


Tp 

^2hk 


rt(k, 


w = — 


W — £’ 2 hk 

Ah(k) 




1/2 

*^ 2 hk 


U + -E2hk ’ 


( 110 a) 


[w-^k-S) ^(k,a;)][a;-h^k + S) '’(k,-w)] - A2(k) 
aikAh(k) 


-h 


2-Elhk 

a2kAh(k) 


1 


1 


w — -Eihk 

1 


^ + -^'Ihk 

1 


2 -E 2 hk Vo; — £' 2 hk ^ -|- ill 2 hk 


( 110 b) 


where aik = {El^^ - Ml)/{El^^ - £ 3 ^^), a 2 k = (^|hk “ ^k)/(^i\k “ and 

there are four branches of the charge-carrier quasiparticle spectrum due to the 
presence of both the normal-state pseudogap and SC gap, -Eihkj —-E'lhkj -E 2 hk) and 
-•E 2 hk, with Eihk = -v/iflk + 0 k]/ 2 , £^ 2 hk = i/fOk - 0 k]/ 2 , and the functions. 


— ?k +-^k + 8Apg(k)-|-A^(k), 

ek = - M2)/3ik + 16A2g(k)/32k + A4(k), 


( 111 a) 

( 111 b) 


with /3ik = ^k “ + 2A^(k), / 32 k = ('^k — Afk)^ + A^(k), while the charge-carrier 

coherence factors. 


Tp 

'/ihk 


T/2 

•^Ibk 


Tp 

'/2hk 


aik 


aik ( 1 - 

-nihk / 


— ask 1 + 


Mk \ 
-E^lhk/ 


ask 


t/2 __ 

'^2hk — 


a2k 


^^^2hk) 


ask 


-fit) 


a2k 1 — 


Ck 

-E2hk 


— ask 1 — 


M,< 


-1 


-E2hk/ 


( 112 a) 

( 112 b) 

( 112 c) 

( 112 d) 


satisfy the sum rule t^ihk+^ihk + ^ 2 hk+^ 2 hk ~ ^ wave vector k (normaliza¬ 

tion), where ask = [ 2 Apg(k)]^/(i?^[^j^ — i3|j^j^), and the corresponding normal-state 
pseudogap Apg(k) and energy spectrum Mk i n Eq. (108) can be obtained explicitly 
in terms of the self-energy S^*^^(k, w) in Eq. (108) as. 


^pg(k) — 


Afk = 


L2(k) 

2 /^’ 

L2(k) 

Li(k)’ 


(113a) 

(113b) 


with the functions Ei(k) = —^[^^(kja; = 0 ) and E 2 (k) = S^^^(k, w = 0 ), and can 

be obtained directly from the charge-carrier self-energy S^*^^(k, w) in Eq. (62a). 
In this case, it is t hen straightforward to obtain the normal-state pseudogap pa¬ 
rameter from Eq. (113) as |92], 


^Pg- AT 




(114) 
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In comparison with the charge- carrie r Green’s functions in Eq. (60), the charge- 
carrier Green’s functions in Eq. (110) are still BGS-like with the d-wave symmetry, 


although the charge-carrier quasiparticle spectrum has been split further due to 
the presence of the normal-state pseudogap. However, it should be noted that in 
spite of the origins of the normal-st ate p seudogap state, the main feature of the 
charge-carrier propagators in E q. (110) is very similar to these proposed from 
the preformed pair theory 212 , where the pseud ogap state is associated with the 
preformed pairs, and then the calculated result 218 of the conductivity in the 
underdoped cuprates is consistent with the experimental observations, or from 

2Ti\ of the normal-state pseudogap state 


several phenomenological theories 219 


based on the d-wave BGS formalism, where by introducing a phenomenological 
doping and temperature dependence of the normal-state pseudogap, the two-gap 
feature in cupra te superconductors is reproduced, or from a phenomenological 
theory [222 , 223 of the normal-state pseudogap state based on the RVB theory, 
where an ansatz is proposed for the coherent part of the single particle Green’s 
function in a doped RVB state, and then the calculated result of the electronic 
properties in the normal-state pseudogap phase is in qualitative agreement with 
the experimental data. 


6.3. Doping and temperature dependence of normal-state pseudogap 
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Figure 24: The normal-state pseudogap parameter (2Apg) (solid line) and charge-carrier 
pair gap parameter (2Ah) (dashed line) as a function of doping with T = 
0.002J for t/J = 2.5, t'/t = 0.3, and J = llOmeV. Inset: the corresponding 
experimental data of cuprate superconductors taken from Ref. [16| . [From 
Ref. 
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The normal-state pseudogap parameter Apg in Eq. (114) has been calculated 


92 217 , and the results of 2Apg (solid line) and charge-carrier pair gap parameter 


2Ah (dashed line) as a function of doping with T = 0.002J for tj J = 2.5, t'/t = 0.3, 
and J = llOmeV are shown in Fig. 24 in comparison with the corresponding 
experimental data of cuprate superconductors (inset). Obviously , the twO' 
gap feature observed on different families of cuprate superconductors 
is qualitatively reproduced. 


Although the charge-carrier pair gap parameter ^ 


has a domelike shape of the doping dependence, the magnitude of the normal-state 
pseudogap parameter Apg is particularly large in the underdoped regime, and then 
smoothly decreases with increasing doping, eventually disappearing together with 
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superconductivity at the end of the SC dome. In particular, since this normal- 
state pseudogap Apg is directly rela ted to the charge-carrier quasiparticle coherent 
weight ZhF as shown in Eq. (109), the decrease of Apg with increasing doping 
leads to that the charge-carrier quasiparticle coherent weight near the charge- 
carrier Fermi surface grows linearly with doping, which together with the charge- 
carrier pair gap parameter Ah shows that only S number of the coherent doped 
carriers are recovered in the SC-stat e, consistent with the picture of a doped Mott 
insulator with S charge carriers 


4,5 53 138 


ilp 


224 . This is much different from the 


case in conventional superconductors [34 ,35 , where the charge-carrier coherent 
weight Zp ~ 1 near the Fermi surface, since the normal-state of conventional 
superconductors is a standard Landau Fermi-liquid. 



Figure 25: The normal-state pseudogap parameter as a function of temperature at d = 
0.09 for tjJ = 2.5 and t'jt = 0.3. [From Ref. 


136 
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Figure 26: The normal-state pseudogap crossover temperature T* (solid line) and super¬ 
conducting transition temperature (dashed line) as a function of doping for 
tj J = 2.5, i' jt — 0.3, and J = llOmeV. Inset: the corresponding experimen¬ 
tal data of cuprate superconductors taken from Ref. 1^. [From Ref. 
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As in the temperature dependence of the charge-carrier pair gap (see Fig. [^, 
the normal-state pseud ogap is also temperature dependent. To show this point 
clearly, the result [136 of the normal-state pseudogap parameter as a function of 


temperature at h = 0.09 for tjJ = 2.5 and t'/f = 0.3 is shown in Fig. 25 where in 


analogy to the temperature dependence of the charge-carrier pair gap parameter 
(see Fig. [^, the normal-state pseudogap parameter decreases with increasing 
temperatures, and vanishes when temperature reaches the normal-state pseudogap 
crossover temperature T *, then the system crossovers from the strongly correlated 
normal-state pseudogap state at the temperatures T < T* to the normal-metal 
phase with largely transport properties at the temperatures T > T*. Furthermore, 
this normal-state pseudogap crossover temperature T* is also doping dependent, 
and the results of T* (solid line) and Tc (dashed line) as a function of doping 


for t/ J = 2.5, t'/t = 0.3, and J = llOmeV are shown in Fig. 26 in comparison with 
the experimental data obtained from different families of cuprate superconductors 


In corresponding to the results of the doping dependence of Apg and Ah 


m 


24 


T* is much larger than Tc in the underdoped regime, and then smoothly 
decreases with increasing doping. Moreover, both T* and Tc converge to the end 


of the SC dome, in qualitative agreement with the experimental results 14-22 
[ 30 ] . In particular, in comparison with the doping dependence of the couphng 
(interaction) strength 14^ in Fig. [^a), it is therefore found 92 that all T*, Apg, 
and Ves show the same doping dependence, i.e., T* ~ Apg ~ I4ff, while such 


a relationship among T*, Apg, and I4ff has been confirmed experimentally on 


cuprate supercond uctor s [l2 
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The results ^2 217 in Fig. 


24 and Fig. 26 show clearly that there are two 


coexisting energy gaps in the whole SC dome: one associated with a direct mea¬ 
surement of the binding energy of the two charge carriers forming a charge-carrier 
pair, while the other with a suppression of the spectral weight of the low-energy ex¬ 
citation spectrum. The charge carriers interact by the exchange of spin excitations 
in both the particle-particle and particle-hole channels, with the coupling (inter¬ 
action) strength I4fj that falls linearly with doping from a strong-coupling case in 
the underdoped regime to a weak-coupling side in the overdoped regime (see Fig. 
1^). On the other hand, the normal-state pseudogap Apg(k) originates from the 
self-energy in the particle-hole channel due to the charge-carrier interaction, which 
leads to that the normal-state pseudogap parameter Apg and normal-state pseu¬ 
dogap crossover temperature T* follow qualitatively a doping dependence in like 
manner of the coupling strength I4ff- In particular, since the charge-carrier inter¬ 
actions in both the particle-hole and particle-particle channels are mediated by the 
same spin excitations as we have shown in Eq. (54), all these charge-carrier interac¬ 
tions in the particle-hole and particle-particle channels are controlled by the same 
magnetic exchange coupling J. In this sense, both the normal-state pseudogap 
and SC gap in the phase diagram of the cuprate superconductors are dominated 
by one energy scale. Moreover, the theory of the kinetic-energy driven supercon¬ 
ductivity starting from the t-J model (34) also shows that both the normal-state 
pseudogap and SC gap in the cuprate superconductors are the result of the strong 
electron correlation, a generic consequence of the strong Coulomb interaction in 
the large-17 Hubbard model (Q, and they have the same magnetic origin. 


6.4. Summary 

In this section, we have given a brief review of the microscopic theory of the normal- 
state pseudogap state in cuprate superconductors. The charge-carrier interaction 
arising through the exchange of spin excitations that generates charge-carrier pair 
state in the particle-particle channel also induces the normal-state pseudogap state 
in the particle-hole channel, therefore the kinetic-energy driven SC mechanism 
provides a natural explanation of both the origin of the normal-state pseudogap 
state and the SC mechanism for superconductivity in cuprate superconductors. 
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In this microscopic theory of the normal-state pseudogap state, the normal-state 
pseudogap has been identified as being a region of the self-energy effect in the 
particle-hole channel in which the normal-state pseudogap suppresses the spectral 
weight of the low-energy excitation spectrum. Furthermore, this microscopic the¬ 
ory of the normal-state pseudogap state also indicates that (a) there is a true order 
parameter defining the existence of the normal-state pseudogap phase; (b) there is 
a coexistence of the SC gap and normal-state pseudogap in the whole SC dome; 
(c) both the normal-state pseudogap and the SC gap are dominated by one energy 
scale, and they are the result of the strong electron correlation; (d) the normal- 
state pseudogap is directly related to the quasiparticle coherence, and therefore 
antagonizes superconductivity; (e) the correct phase diagram with respect to the 
normal-state pseudogap line is that the normal-state pseudogap is much larger 
than that of the SC gap in the underdoped regime, and then it merges gradually 
with the SC gap in the overdoped regime, eventually disappearing together with 
superconductivity at the end of the SC dome. 


7. Charge transport 


Since the discovery of superconductivity in cuprte superconductors a sig¬ 

nificant body of reliable and reproducible data has been accumulate by using 
many probes [6 14-22 


which shows that the most remarkable expression of the 
nonconventional physics is found in the normal-state. The normal-state properties 
in the underdoped and optimally doped regimes exhibit a number of anomalous 
properties in sense that they do not fit in with the standard Landau Fermi-liquid 
theory. However, there is mounting evidence that the anomalous normal-state 
properties in the under doped an d optimally doped regimes is dominated by the 

Among the striking features of the normal-state 


22 


normal-state pseudogap |6|14 

properties in the underdop^ and optimally doped regimes, the physical quantity 
which most evidently displays the signature for the normal-state pseudogap is the 


charge transport 14,76 


225-235 , which is manifested by the conductivity and 


resistivity. The opticaTstudies of the quasiparticle excitations have revealed much 
about the nature of the charge carriers in cuprate superconductors. In particular, 
the normal-state pseudogap can be observed directly by the infrared measurements 
of the conductivity. Experimentally, it has been shown in terms of the Kramers- 
Kronig analysis of the reflectance that th e conduc tivity is rather universal within 
the whole cuprate superconductors [l4|76p% - 235| , where a key feature is the two- 
component conductivity: a narrow band centered around energy cj ~ 0 followed by 
a broadband centered in the midinfrared region in the underdoped and optimally 
doped regimes. The conductivity in the underdoped and optimally doped regimes 
shows a non-Drude behavior (the conductivity decays as —)■ l/cu) at low energies, 
and is carried by 6 charge carriers, while the midinfrared spectral weight is biased 
towards the low-energy region with the increase of doping. In particular, this two- 
component conductiv ity exten ds to the normal-state pseudogap boundary in the 
phase diagram at T* [231 - 234 . 

Within the framework of the kinetic-energy driven SC mechanism, the doping 
and temperature dependence of the conductivity have been studied in the whole 
doping range from the underdoped to heavil y overdoped by considering the effect 
of the normal-state pseudogap 88,236 - 239 , and the result shows that the part 


of the low-energy spectral weight of the conductivity spectrum in the underdoped 
and optimally doped regimes is transferred to the higher energy region to form the 
unusual midinfrared band, however, the onset of the region to which the spectral 
weight is transferred, is always close to the normal-state pseudogap. 
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7.1. Linear response theory 


Through the standard linear response theory 131 , t he hnite- frequency conductiv¬ 
ity of cuprate superconductors can be expressed as (M6H23^ , 


a(uj) = - 


Imn(a;) 


a; 


(115) 


with the electron current-current correlation function, 

n(r-T') =-(T.rj(T) •j(r')), (116) 

where the electron current operator j is obtained by evaluating the time-derivative 
of the polarization operator (70), and has been given explicitly in Eqs. (74) and 
(75) as, 

iex2 t' 
h 


J = 


iexit 

h 




t 

l+f]cr' 


hla + 

lf]ij Ira 




(117) 


In the normal-state, the electron current-current correlation function is evaluated 
in terms of the full charge-carrier Green’s function as. 


with the current vertex, 


1 , 


7sk = -^[{xit - ‘^X 2 t' cos kyf s\v? kx + {xit - cos kxf SIT? ky], (119) 

while the full charge-carrier Green’s function in the normal -state can be obtained 
in terms of the full charge-carrier diagonal Green’s function (110a) in the con ditio n 


of the charge-carrier pair gap Ah = 0, and has been evaluated explicitly as 238 

+ 


a 


(n) 

2k 




UJ - E, 


( 120 ) 


where there are two branches of the charge-carrier quasiparticle spectrum due to 
the presence of the normal-state pseudogap. 


F+ - ^ 

-^hk “ o 


^hk- 2 


Ck - Mk + J(ek + Mk)2 + 16A2 (k) 


4 - Mk - J(4 + Mk)2 + 16A2 (k) 


(121a) 

(121b) 


while = {E^k +/{E+^-E^^J and = -(^hk + ^k)/(£^hk--E^hk) satisfy 

the sum rule = 1 for any wave vector k. In this case, the conductivity 

in the normal-sate is obtained in terms of the full charge-carrier Green’s function 
(pOl) as. 




7sk 


do;' , ,, , s , ,.np(u}' + uj) — nF(uj') 

—7lh(k,a;' + a;)7lh(k,a;')^^- ^(122) 


27r 


OJ 


with the charge-carrier spectral function Ah(k, w) = —2Img(k,a;). In the fermion- 
spin theory, although both char ge car riers and spins contribute to the charge dy¬ 
namics, the results in Eqs. (117) and ( 122[ ) sho w that th e anomalous conductivity 
properties are mainly caused by charge carriers |236 - 239 , which are strongly renor¬ 
malized because of the strong interaction with the fluctuation of the surrounding 
spin excitations, and then the low-energy spectral weight of the conductivity spec¬ 


trum is proportional to the charge-carrier doping concentration 5 226 227 
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7.2. Doping dependence of conductivity in normal-state 



Figure 27: The conductivity in the normal-sate as a function of energy a.t 6 = 0.09 (solid 
line), S = 0.15 (dashed line), and S — 0.25 (dotted line) with T = 0.002J for 
t/J = 2.5 and f /t = 0.3. [From Ref. 


238 


T he condu ctivity < 7 ( 0 ;) in Eg. (122) in the normal-state has been calcu¬ 
lated [236 -238 , and the result [238 of < 7 ( 0 ;) as a function of energy at <5 = 0.09 
(solid line), 6 = 0.15 (dashed line), and 6 = 0.25 (dotted line) with T = 0.002J 
for t/J = 2.5 and t'/t = 0.3 is shown in Fig. \Tl\ hereafter in this section we 
set the charge e and reduced Planck constant h as the unity. Obviously, this cal¬ 
culated result captures all qualitative features of the doping dependence of the 
conducti vity observed experimentally on cuprate superconductors in the normal- 
state [14 


76,225-235 . In the underdoped regime, there are two bands in the 


conductivity aioj) separated by a gap at a; ~ 0.2t. The higher energy band, corre¬ 
sponding to the midinfrared band, shows a broad peak at w ~ 0.38t. In particular, 
the transferred weight of the low-energy band forms a sharp peak at ui ~ 0, which 
can be described formally by the non-Drude formula, while the onset of the re¬ 
gion to which the spectral weight is transferred, is always close to the normal-state 
pseudogap Apg, reflecting a fact that due to the presence of the normal-state pseu¬ 
dogap in cuprate superconductors, the part of the low-energy spectral weight of the 
conductivity spectrum in the normal-state in the underdoped regime is transferred 
to the higher energy region to form the unusual midinfrared band. In the other 
words, the appearance of the higher energy midinfrared band is closely related 
to the effect of t he normal-state pseud ogap on the infrared response in cuprate 

In a given doping, the spectral weight is 

However, the 


231i|233i|234[|240[|241 

ne area under the conductivity curve in Fig. 


27 


superconductors 
proportional to t 

weight and position of the midinfrared band are strongly doping dependent. In 
particular, the result in Fig. indicates that as the charge-carrier doping in¬ 
creases, although the overall conductivity increases, the magnitude of the gap in 
the conductivity spectrum decreases, and then the higher energy midinfrared band 
moves towards the low-energy non-Drude band. In the optimal doping, although 
two band features are still apparent, the positions of the gap in the conductivity 
spectrum and midinfrared peak appreciably shift towards the lower energies at 
oj ~ 0.16t and uo ~ 0.3t, respectively, reflecting a tendency that with increasing 
doping, the magnitude of the gap in the conductivity spectrum decreases, while 
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and 
AJu„ with 


the midinfrared band moves towards the low-energy non-Drude band. However, as 
in the case in the underdoped regime, the low-energy peak in the optimal doping 
still shows the non-Drude formula. This follows from a fact that the re sult o f the 
conductivity spectrum in the optimal doping 5 = 0.15 has been fitted 238| 
the fitted result shows that the lower-energy peak decay perfectly as 
A ~ 0.01. On the other hand, the tendency of the decrease of the magnitude of the 
gap in the conductivity spectrum and the midinfrared band moving towards to the 
low-energy non-Drude band with increasing doping is particularly obvious in the 
overdoped regime. In particular, the low-energy non-Drude peak incorporates with 
the midinfrared band in the heavily overdoped regime, and then the midinfrared 
feature disappears, which leads to that the low-energy Drude type behavior of the 
conductivity recovers, and then in contrast to t he c ase in the underdoped and 


238 


optimally doped regimes, the lower-energy peak 
the heavily overdoped regime, with A ~ 0.01 and B 


decay as —)■ H/-|- B) in 
' 0.0045. In sectionit has 
been shown that the magnitude of the normal-state pseudogap (then T*) is partic¬ 
ularly large in the underdoped regime, and then smoothly decreases upon increase 
of doping. The calculated result of the doping dependence of the conductivity 
in Fig. 27 also implies that the onset of the region to which the spectral weight 


is transferred (then the midinfrared peak) shows the same trend with doping in 
like manner of the doping dependence of the normal-state pseudogap (then T*). 
To show this point clearly, we plot the position of the midinfrared peak (dashed 
line) and T* (solid line) as a function of doping with T = 0.002J for t/ J = 2.5 
and t '/t = 0.3 in Fig. 28 in comparison with the corresponding experimental 
data [M] of YBa 2 Cu 30 y(inset). This calculated result in Fig. is very well 


consistent with the experimental data observed from the conductivity measure¬ 
ments [M ,231 , and therefore confirm the effect of the normal-state pseudogap on 
the infrared response. 





* 



5 


Figure 28: The position of the midinfrared peak (dashed line) and T* (solid line) as a 
function of doping with T = 0.002J for t/J = 2.5 and t' jt = 0.3. Inset: the 
corresponding experimental data of YBa 2 Cu 30 y taken from Ref. 


231 


The low-energy non-Drude peak and unusual midinfrared band of the con¬ 
ductivity spectrum in the normal-state in the underdope d an d optimally doped 
regimes are also temperature dependence. In Fig. 29 we [238] show the conduc¬ 
tivity in the normal-state as a function of energy at 0 = 0.09 with T = 0.02 J (solid 
line), T = 0.146J (dashed line), and T = 0.186J (dotted line) for t/J = 2.5 and 
t'/t = 0.3. In section]^ the calculated normal-state pseudogap crossover temper¬ 
ature is T* ~ 0.19J at <5 = 0.09. The result in Fig. [2^ shows that the weight 
of the midinfrared band is severely suppressed with increasing temperatures. In 
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Figure 29: The conductivity in the normal-sate as a function of energy at S = 0.09 with 
T = 0.02J (solid line), T = 0.146J (dashed line), and T = 0.186J (dotted 
line) for tjJ = 2.5 and t'/t = 0.3. [From Ref. 


238 


particular, the weight of the midinfrared band vanishes above the temperatures 
T > T*, and then the low-energy Drude type behavior of the conductivity is re¬ 
covered, which is also qualitatively consistent with the experimen tal data observed 
in cuprate superconductors in the normal-state 76 225-230 . 


7.3. Doping dependence of conductivity in superconducting-state 

The discussions of the conductivity of cuprate superconductors in the normal-state 


in subsection 7.2 can also be generalized to the SC-state in the condition of the 
charge-carrier pair gap Ah ^ 0. In the SC-state, the electron current-current 
correlation function can be obtain ed in ter ms o f the full charge-carrier diagonal 
and off-diagonal Green’s functions (110) as |239|, 


U{iujn) = -^(^e)2^^7sk^ +ia;„)5r(k,ia;n/) 

k 

r(k,za;„/-I-fa;„)r'l'(k,iwn')]) 


(123) 


and then the finite-frequency conductivity (|115 ) of cuprate superconductors in the 


SC-state can be evaluated explicitly as 239 
2 


1 _ r°° A,,! 

~J[J ^[^ 3 (k,w + a;')7lg(k,a;') + Ar(k,w + a;')7lr(k,a;')] 


np^uj') — np{uj + uj') 


UJ 


(124) 


where the spectral functions Ag(k,t(j) and Ar(k,cj) are obtaine d in terms of the 
charge-carrier diagonal and off-diagonal Green’s functions in Eq. (110) as Ag(k, oj) = 
— 2 Im/ 7 (k,a;) and ^ r(k, c u) = —2Imr^^(k, w), respectively. 


In Fig. 30 


we 


239] show the calculated result of the conductivity (124) 


m 


the SC-state as a function of energy at d = 0.09 for t/J = 2.5 and t'/t = 0.3 
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Figure 30: The conductivity in the SC-state as a function of energy at <5 = 0.09 with T = 
0.002J for t/J = 2.5 and t'/t = 0.3. Inset: the corresponding experimental 
data of the underdoped Bi2Sr2CaCu208+5 taken from Ref. |233| . [From Ref. 


with T = 0.002J in comparison with the corresponding experimental result [233 


of the underdoped Bi2Sr2CaCu208+5 (inset). The result in Fig. [30| shows clearly 
that the two-component feature of the conductivity spectrum iiiThe SC-state is 
the same as in the normal-state case. In particular, in comparison with the result 
of the conductivity in the normal-state in subsection 7.2, it is found that the 


spectral weight of the low-energy component of the conductivity in the SC-state 
is further suppressed by the SC gap, however, there is no depletion of the spectral 
weight of the higher energy midinfrared band of the conductivity in the SC-state, 
whic h is consistent with the experimental observation on cuprate superconductors 
[231 . In this case, although there is a coexistence of the SC gap and normal-state 


pseudogap, the onset of the region to which the spectral weight is transferred is 
also close to the normal-state pseudogap Apg, then in analogy to the evolution 
of the conductivity spectrum with doping in the normal-state, the positions of 
the gap in the conductivity spectrum a nd m idinfrared peak gradually shif t to the 
lower energies with increasing doping 239j. All the calculated results [239 of 


the conductivity spectrum in the SC-state are also qualitatively consistent with 
the corresponding experimental data of cuprate superconductors in the SC-state 
[l4l[76l [22$j^ . 


In the SC-state, there are two parts of the charge-carrier quasiparticle contri¬ 
bution to the redistribution of the spectral weight in the conductivity spectrum 
in the SC- state: the contribution from the first term of the right-hand side in 
Eq. (124) comes from the spect ral fu nction obtained in terms of the charge- 
carrier olagonal Green’s function (llOa), and therefore is closely associated with 
the normal-state pseudogap Apg in the particle-hole channel, whil e the additional 
contribution from the second term of the right-hand side in Eq. (124) originates 
from the spectral funct ion obtained in terms of the charge-carrier off-diagonal 
Green’s function (110b), and is closely related to the charge-carrier pair gap Ah 
in the particle-particle channel. However, since Ah Apg in the underdoped and 
optimally doped regimes as we have mentioned in section the charge-carrier pair 
gap only suppresses the spectral weight of the low-energy component, while the 
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normal-state pseudogap related shift of the spectral weight from the low-energy to 
the higher energy midinfrared band in the SC-state conductivity spectrum becomes 
arrested. 

The effect of the normal-state pseudogap on the infrared response in cu prate 
superconductors has been also discussed based on the preformed pair th eory 
and the phenomenological theory of the normal-state pseudogap state 


242 


the results of the unusual two-component conductivity spectrum are qualitatively 
consistent with the above obtained result based on the kinetic-energy driven SC 
mechanism by considering the interpla y betwee n the SC gap and normal-state 
pseudogap. In particular, their results 218,242 also indicate that in the under¬ 


doped and optimally doped regimes, the transfer of the part of the low-energy 
spectral weight of the conductivity spectrum to the higher energy region to form a 
midinfrared band is intrinsically associated with the presence of the normal-state 
pseudogap. 

In a standard Landau Fermi-liquid, the shape of the conductivity spectrum 
cr{uj) is normally well accounted for by the low-energy Drude formula that de¬ 
scribes the free charge carrier contribution to cr{uj), and then when the tempera¬ 
tures T < Tc, the spectral weight of the condensate in the SC-state comes from 
low energies 1^. However, in cuprate superconductors, the part of the low-energy 
spectral weigHt in the conductivity spectrum in the underdoped and optimally 
doped regimes is transferred to the higher energy region to form the unusual mid¬ 
infrared band, and then the width of the low-energy band is narrowing, while 
the onset of the region to which the spectral weight is transferred is close to the 
normal-state pseudogap Apg. Moreover, since the unusual midinfrared band is 
taken from the low-energy band, so that both the low-energy non-Drude peak (the 
conductivity decays as —)■ l/w at low energies) and unusual midinfrared band de¬ 
scribe the actual charge-carrier density. In the framework of the kinetic-energy 
driven SC mechanism, the normal-state pseudogap state is the result of the strong 
electron correlation, and therefore the transfer of the part of the low-energy spec¬ 
tral weight of the conductivity spectrum in the underdoped and optimally doped 
regimes to the higher energy region to form the unusual midinfrared band is a 
natural consequence of the strongly correlated nature in cuprate superconductors. 
In particular, this strong electron correlation which induces a shift of the spectral 
weight from the low-energy to the higher energy midinfrared band in the con¬ 
ductivity spectrum, has been confirmed by the early numerical simulations based 

■S and in the SC-state |85l. In 


on the t-J model in the normal-state 101 243 


the normal-state, the 1/co perfect decay of the conductivity at low energies in the 
optimally doped regime is closely related with the linear temperature resistivity, 
since it reflects an anomalous frequency dependent scattering rate proportional to 
CO instead of co'^ as would be expected in the standard Landau Fermi-liquid. This 
linear temperature resistivity is one of the characteristically anomalous properties 
of cuprate superconductors in the normal-state, and ha s bee n also phenomenolog¬ 
ical discussed within the marginal Fermi-liquid theory 246 . In particular, based 


on the slave-boson gauge theory, it has been shown within the t-J model that 
above the Bose-Einstein temperature, the boson inverse lifetime due to scattering 
by the gauge field is of order T, which suppresses the condensation temperature 
and leads to a linear T resistivity ro6,67 . However, in the SC-state, the large 


normal-state pseudogap in the underdoped and optimally doped regimes heavily 
reduces the fraction of the charge carriers that condense in the SC-state 


7.4. Summary 

The calculated result of the conductivity spectrum summarized in this section 
shows very clearly that if the effect of the normal-state pseudogap is taken into 
account in the framework of the kinetic-energy driven SC mechanism, the conduc¬ 
tivity of the t-J model calculated based on the linear response approach per se can 
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correctly reproduce the main features found in infrared response measurements 
on cuprate superconductor in both the normal- and SC-states. The conductivity 
spectrum in the underdoped and optimally doped regimes contains the low-energy 
non-Drude peak and unusual midinfrared band. However, the position of the mid- 
infrared band shifts towards to the low-energy non-Drude peak with increasing 
doping. In particular, the low-energy non-Drude peak incorporates with the mid- 
infrared band in the heavily overdoped regime, and then the low-energy Drude 
behavior recovers. The qualitative reproduction of the main features of infrared 
response measurements on cuprate superconductors also shows that the transfer 
of the part of the low-energy spectral weight in the conductivity spectrum in the 
underdoped and optimally doped regimes to the higher energy region to form 
the unusual midinfrared band can be attributed to the effect of the normal-state 
pseudogap on the infrared response in cuprate superconductors. 


8. Conclusion and discussions 


In this article, we have given a brief review of the kinetic-energy driven SC mech- 
where the main conclusions are summarized as: 


anism. 


(a) In the fermion-spin theory (23), the constrained electron is decoupled as a 
product of a charge carrier and a localized spin, and the charge carrier represents 
the charge degree of freedom of the constrained electron together with some effects 
of the spin configuration rearrangements due to the presence of the doped charge 
carrier itself, while the spin operator represents the spin degree of freedom of the 
constrained electron. In the decoupling scheme, this fermion-spin representation 

a is a natural representation of the constrained electron dehned in a restricted 
ert space without double electron occupancy. The main advantage of the 
fermion-spin theory ( |23[ ) is that the electron local constraint for single occupancy 
is satisfied in actual calculations. In particular, the charge carrier or spin itself 
is U{1) gauge invariant, and in this sense, the collective modes for the charge 
carrier and spin are real and can be interpreted as the physical excitations of 
cuprate superconductors. Although both charge carriers and spins contribute to 
the charge and spin dynamics, the charge-carrier relaxation time is responsible to 
the charge transport, and the spin relaxation time is responsible to the dynamical 
spin response, while as a result of the charge-spin recombination, the the electronic 
properties are dominated by electron quasiparticles. This is an efficient calculation 
scheme which can provide very good results even at the MF level. 

(b) In the framework of the kinetic-energy driven SC mechanism developed 
based on the fermion-spin theory (23), the charge-carrier pairing interaction orig¬ 


inates directly from the kinetic energy of the t-J model (34) by the exchange of 


spin excitations in the higher powers of the doping concentration, and then these 
charge-carrier pairs (then the electron Cooper pairs) condense to the d-wave SC- 
state, where the spin excitation in cuprate superconductors plays a similar role to 
that of the phonon in conventional superconductors. Although the physical prop¬ 
erties of cuprate superconductors in the normal-state are fundamentally different 
from these in the standard Landau Fermi-liquid state, the kinetic-energy driven 
SC-state still is conventional BCS-like with the d-wave symmetry, and then the 
obtained formalism for the charge-carrier pairing can be used to compute Tc and 
the related SC coherence of the low-energy excitations in cuprate superconductors 
on the hrst-principles basis much as can be done for conventional superconductors. 
Moreover, the kinetic-energy driven charge-carrier pair state is controlled by both 
the charge-carrier pair gap and quasiparticle coherence, which leads to that the 
maximal Tc occurs around the optimal doping, and then decreases in both the un¬ 
derdoped and overdoped regimes. This kinetic-energy driven SC mechanism also 
indicates that the strong electron correlation favors superconductivity, since the 
main ingredient is identified into a charge-carrier pairing mechanism not from the 
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external degree of freedom such as the phonon, but rather solely from the internal 
spin degree of freedom of the constrained electron. 

(c) The same charge-carrier interaction arising through the exchange of spin 
excitations that induces the d-wave SC-state in the particle-particle channel also 
generates the normal-state pseudogap state in the particle-hole channel, therefore 
there is a coexistence of the SC gap and normal-state pseudogap in the whole SC 
dome. Consequently, this normal-state pseudogap is identified as being a region of 
the self-energy effect in the particle-hole channel in which the normal-state pseu¬ 
dogap suppresses the spectral weight of the low-energy excitation spectrum. This 
normal-state pseudogap vanishes at the normal-state pseudogap crossover temper¬ 
ature T*, with T* that is much larger than Tc in the underdoped and optimally 
doped regimes, and monotonically decreases upon the increase of doping, even¬ 
tually disappearing together with Tc at the end of the SC dome. In particular, 
the normal-state pseudogap is directly related to the quasiparticle coherence, and 
therefore antagonizes superconductivity. Moreover, both the normal-state pseu¬ 
dogap and charge-carrier pair gap are dominated by one energy scale, and they 
are the result of the strong electron correlation in cuprate superconductors, while 
the domelike shape of the doping dependence of Tc, the monotonic decrease of T* 
with doping, and relatively anomalous normal-state properties are a natural con¬ 
sequence of the Mott physics in which double occupancy is suppressed by strongly 
Coulombic repulsion. The theory also indicates that the kinetic-energy driven SC 
mechanism provides a natural explanation of both the origin of the normal-state 
pseudogap state and the charge-carrier pairing mechanism for superconductivity. 

(d) Within the framework of the kinetic-energy driven SC mechanism, a num¬ 
ber of typical properties of cuprate superconductors have been studied. In this 
review article, the selected results are summarized, including the doping depen¬ 
dence of the electromagnetic response, the dynamical spin response from low- 
energy to high-energy, and the charge transport, and are qualitatively comparable 
to the corresponding experimental results observed in cuprate superconductors. 
Furthermore, this kinetic-energy drive n BC S-type formalism gives an explana¬ 
tion of the Raman scattering spectra |247| obtained in terms of the electronic 


Raman response measurement techni que 
tion of the thermodynam ic properties 
ments 


162 163 


248-251 . In particular, it 


number of the sC-state properties in the presence of impurities 252 


?9|. It also gives a consistent descrip- 
rserved from heat capacity measure- 
las been used to suc cessf ully describe a 

including the 


microwave conductivity 253 , the scanning tunneling microscopic measurements 


of the coherent Bogoliubov quasiparticle dispersion, and the related ex tinct ion of 
Bogoliubov quasiparticle scattering interference at low temperatures [254 . Es¬ 
tablishing these agreements between the calculated results obtained based on the 
kinetic-energy driven SC mechanism and the corresponding experimental data ob¬ 
served from a wide variety of measurement techniques are important to confirm 
the nature of the SC phase of cuprate superconductors to be the kinetic-energy 
driven d-wave SC-state. 

In this review article, we have restricted our attention to the hole-doped cuprate 
superconductors. However, supercon ductivity in cuprates also emerges when elec¬ 
trons are doped into Mott insulators [255 ,256 . Both the hole-doped and electron- 
doped cuprate superconductors have the layered struct ure of th e square lattice 
of the Cu02 plane separated by insulating layers 1-31255 256 . Although the 


signihcantly different behaviors of the hole-doped aim dectron-doped cuprate su¬ 
perconductors are observed due to the electron-hole a symm etry, the symmetry of 
the SC order parameter is common in both case 39,257 , manifesting that two 
systems have similar underlying SC mechanism. In particular, the strong elec¬ 
tron correlation is common for both the hole-doped and electron-doped cuprate 
superconductors, and then it is possible that superconductivity in the electron- 
doped cuprate superconductors is also driven by the k inetic energy as in the hole- 
doped case. In this case, the charge asymmetry [258 in superconductivity of the 
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electron- doped and hole-doped cuprate superconductors and the electronic Raman 
response 259| in the electron-doped cuprate superconductors have been discussed 
based on the kinetic-energy driven SC mechanism, and the calculated results are in 
qualitative agreement with the experimental data observed on the electron-doped 
cuprate superconductors. 

Besides the square lattice cuprate superconductors, some cuprate materials 


262 , such as Sri 4 Cu 2404 i, do not contain Cu 02 planes common to cuprate 


superconductors but consist of two-leg CU 2 O 3 ladders and edge-sharing Cu 02 
chains. In particular, the doped two-leg ladder cuprates are a system in which 
the SC state i s realized by applying a high pressure in the highly charge-carrier 
doped regime 261, 262|. These ladder cuprate materials also are natural exten¬ 


sions of the Cu-U chain compounds towards the Cu 02 sheet structures. Within 
the framework of the kinetic-energy driven SC mechanism, some typical proper¬ 
ties of the two-leg ladder cupr ate supe rconductors have been stud ied, including 
the pressu re depend ence of Tc [263, 264 , the charge dynamics [265 , and the spin 


dynamics [266 267 , and the calculated results are also qualitatively consistent 


with the experimental data obtained from the experimental measurements on the 
two-leg ladder cuprate superconductors. 

Finally, we should be noted that much remains to be done. In particular, for the 
normal-state pseudogap, which grows upon underdoping, it seems natural to seek 
a connection to the physics of the AF insulating parent compounds 14-^. How¬ 
ever, at half-filling, the t-J model is reduced as the AF Heisenberg rnodel with 
an AFLRO. Although a small density of charge carriers is sufficient to destroy 
AFLRO, this AFLRO remains until the extremely low-doped regime (5 < 0.045) 

(zmiiiii 


As we have mentioned in Eq 
at the extremely low-doped regime, anc 


(|3^, the conduct is disrupted by AFLRO 
then an important issue is how to ex¬ 


tend the theory of the normal-state pseudogap state for the doped regime without 
AFLRO to the case at the extremely low-doped regime with AFLRO for a proper 
description of the connection between the finite doping normal-state pseudogap 
and the zero-doping quasiparticle dispersion. 
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Appendix A. Matrix representation of projection operator 


of the charge-carrier states are given by. 


The charge-carrier operators and hi„ in the basis 


87 






0 0 
1 0 


(Al) 


(A2) 
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while the spin raising and lowering operators S'j^ and in the spin 1/2 space, 


are given by, 







s; 


0 0 \ 

1 oj • 


In the product space [charge) (8) [spin), the basis vectors are (^, 


|l,t) = 


| 14 ) = 


|o,t) = 


|o,i) = 



(A3) 


(A4) 


(A5a) 


(A5b) 


(A5c) 


(A5d) 


which form a complete set, then the fermion-spin transformation defined by Eq. 
(23) in this basis gives the following matrix representation for the constrained 
electron operators [87], 




'•IfSr 



0 0 
1 0 


Cii 







0 0 \ 

1 oA 


0 1 
0 0 


cl = 


0 0 \ 

1 oA 


0 0 
1 0 


/ 0 0 
0 0 
0 0 

V 0 0 

/ 0 0 

0 0 
0 0 

V 0 0 

/ 0 0 

0 0 
0 1 

V 0 0 

/ 0 0 

0 0 
0 0 
Vl 0 


0 0 \ 
1 0 
0 0 
0 0 / 

0 1 \ 
0 0 
0 0 
0 0 / 

0 0 \ 
0 0 
0 0 
0 0 / 

0 0 \ 
0 0 
0 0 
0 0 / 


(A6a) 


(A6b) 


(A6c) 


(A6d) 
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However, as we have mentioned in Eq. the restricted Hilbert space without 
double electron occupancy in the t-J model Q consists of three states, |0), | t); 
I 1), namely. 


| 0 ) = 





(A7) 


To remove the extra degrees of freedom in the |charge) ® |spin) space, we 
introduce a projection operator P. By requiring H|l,t) = -Pllji) = |0), -P|0,t 
I t)) and -P|0,1) = I I), we can easily obtain its matrix representation, 



P = {P.a] 


and its hermitian conjugation. 


= 


/1 

1 

0 

= 0 

0 

1 

V 0 

0 

0 


/1 

0 

0 \ 

1 

0 

0 

0 

1 

0 

V 0 

0 

1 / 


(A8) 


(A9) 


Using this projection operator, the electron operators in the restricted Hilbert 
space without double electron occupancy are given by (87|, 


C„ = Pih'„S^P} = f 0 

Cll = Pih^SlP} = f 0 


1 0 \ 

0 0, (?;, = P,h, 

0 0 / 

0 1 \ 

0 0 , cl = PA 
0 0 / 


/ 0 0 0 
5+P/ =10 0 
\ 0 0 0 

/ 0 0 0 
S,-f/ =000 
\ 1 0 0 


(AlOa) 

(AlOb) 


as quoted in Eq. (31). It is then straightforward to verify the operate 
quoted in Eq. (32). In particular, the charge-carrier number operator. 


relations 



1 0 0 \ 

0 0 0 =|0)(0| = -P(|lt)(ltl + IU)(UI)7’^. (All) 

0 0 0 / 


The physical meaning of Eq. (All) is transparent: the empty state should be 
counted only once, not twice. Since the ME treatment of the constraint on average 


doping concentration 5 is imposed directly on — h\hi, the sum rule for the 

constrained electron is satisfied. 
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